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ABSTRACT 

We report the detection of eclipses in LSPM J1112-I-7626, which we find to be a moderately bright 
{Ic = 12.14 ± 0.05) very low-mass binary system with an orbital period of 41.03236 ± 0.00002 days, 
and component masses Mi = 0.395 ± 0.002 M© and M2 = 0.275 ± 0.001 M© in an eccentric (e = 
0.239 ±0.002) orbit. A 65 day out of eclipse modulation of approximately 2% peak-to-peak amplitude 
is seen in /-band, which is probably due to rotational modulation of photospheric spots on one of the 
binary components. This paper presents the discovery and characterization of the object, including 
radial velocities sufficient to determine both component masses to better than 1% precision, and a 
photometric solution. We find that the sum of the component radii, which is much better-determined 
than the individual radii, is inflated by 3.81q5% compared to the theoretical model predictions, 
depending on the age and metallicity assumed. These results demonstrate that the difficulties in 
reproducing observed M-dwarf eclipsing binary radii with theoretical models are not confined to 
systems with very short orbital periods. This object promises to be a fruitful testing ground for the 
hypothesized link between inflated radii in M-dwarfs and activity. 
Subject headings: binaries: eclipsing - stars: low-mass, brown dwarfs 
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1. INTRODUCTION 

Detached, double-lined eclipsing binaries (EBs) pro- 
vide a largely model-independent means to precisely and 
accurately measure fundamental stellar properties, par- 
ticularly masses and radii. In the best-observed systems 
the precision of these can be at the < 1 per cent level, 
and thus place stringent con straints on stellar evolution 
models (e.g. Andersen 1991: iTorres et al.ll20ldl ). 

Despite the ubiquity of M-dwarfs in the solar neigh- 
borhood, their fundamental properties are poorly un- 
derstood. This is particularly the case below 0.35 Mq, 
the boundary at which field stars are thought to become 
fully convective (e.g. iChabrier &: Baraffg I1997D . Fur- 
thermore, observations of many of the best-characterized 
M-dwarf EBs indicate significant discrepancies with the 
stellar models. 

The components of CM Dra, the best characterized 
double-lined system in the fully-convective mass range 
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(lEggen fc S andage"1967l: lLacvl[l977l: IMetcalfe et all[l996l: 
iMoraleset al.. ,2009i) , have radii approximately 5 — 7% 
larger than predicted by theoretical stellar evolution 
models. The observation of inflated radii for M-dwarf 
binaries at the 5 — 10% level, and effective tempera- 
tures 3 — 5% lower than the models predict, is nearly 
ubiquitous among the best-characterize d objects (e.g. 
Torres et al.ll201C , and references therein; [Morales et al.l 



201c jKraus et al.ll20 1?) , although it is interesting to note 
that [Carter et al.l (12011) flnd a smaller inflation for KOI- 
126B and C, members of a fascinating system containing 
a short-period M-dwarf binary orbiting a K-star, where 
the system undergoes mutual eclipses allowing a precise 
determination of the radii of both M-dwarfs. 

Two main hypotheses have been advanced to ex- 
plain the inflated ra dii in low - mass eclipsing bina- 
ries: metallicity (.Berger et al.l l2006( ). and the ef- 
fect of m agnetic activity (the "activ i ty hypothe- 
sis", e.g. iLqpez-Morales fc RibasI [200l iRibasI 120061: 
IChabrier et al.lboOTD . 
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As discussed by iLopez-MoralesI ()2007[ ). one of the main 
difficulties witli metallicity as an explanation is that cur- 
rent models of M-dwarfs yield radii differing only by a 
very small amount (an increase of approximately 3%) 
between [M/H] = —0.5 to 0, the range spanned by the 
available Baraffe et al.l (|1998[ ) models. Unless the EB 
sample contains many objects of very high metallicity, 
this effect does not seem large enoughQ 

Observationally, the difficulty o f determining rnetal- 
licities for M -dwarfs (e. g. iJB onfils ct^ 120051 : 
iWoolf fc Walle rstcin 2006; B ean et al. 2006) renders the 
metallicity hypothesis difficult to test in practice, and 
this is further exacerbated by the double lines and 
rapid rotation in most M-dwarf eclipsing binary spec- 
tra0 Nonetheless, it is inte resting to n ote that CM Dra is 
thought to be metal poor (jViti et al.l 1997, 2002), which 
would increase the size of the radius discrepancy for this 
object, rather than decreasing it, and more generally, 
this is expected to also be true of other eclipsing bina- 
ries in the sample, assuming t hey follow the metallicity 
distribution of field stars (e.g. iCarnev et al.| [r9891. 

Given these lines of reasoning, it seems unlikely that 
metallicity is the only explanation for the infiated radii 
(although it probably plays a role, particularly in cre- 
ating scatter in the effective temperatures). This has 
led to the "activity hypothesis" becoming the preferred 
explanation in the literature for the infiation. It is mo- 
tivated by noting that nearly all EB systems have short 
orbital periods. For example, below 0.35 M0, the longest 
period system is IRX S J154727.5-I-450803 at 3.55 days 
(jHartman et al.ll2011|) . At such short periods, the com- 
ponents are expected to have been tidally synchronized 
and the orbits circularized by field ages. The effect of 
the companion and tidal locking is likely to significantly 
increase the activity levels. The available observational 
evidence supports this, with the eclipsing binaries be- 
low 0.35 Mq all showing signs of high activity and tidal 
effects including synchronous, rapid rotation, large am- 
plitude out of eclipse modulations. Ha emission. X-ray 
emission, and in most cases, circular orbits. 

Magnetic activity cou ld have an effect either 
by inhibiting c onvec tion (jMullan fc MacDonaldl 120011 : 
iChabrier et al.l I2007D . or due to reduced heat flux 
resulting frorn cool photospheric spots. Of these, 
IChabrier et al.l (|2007f ) find the structure of objects be- 
low 0.35 M0 is almost independent of variations in the 
mixing length parameter, which they use to implement 
inhibition of convection, so the latter possibility (inhibi- 
tion of radiation due to the effect of spots) seems more 
likely the explanation in the mass domain of particular 
interest here. It is important to note that spots also in- 
fluence the solution of eclipsing binary light curves. The 
possibility of s uch systematic error s in the observations 
is discussed bv lMorales et al.l (|2010t ). and we shall return 
to it later in the analysis. 

An obvious way to test the activity hypothesis is to 
measure radii for components of eclipsing binaries with 
long orbital periods, where tides are unimportant, and 

^ As discussed by IL6pez-Morale3 1)20071 ) . it is possible the effect 
on the radius is suppressed by a missing source of opacity in the 
models. Such effects are known to exist as the same models fail 
to reproduce the ob served optical spectra and colors of M-dwarfs 
IIBaraffe et al.|[T99g) . 

^ We thank the referee for bringing this point to our attention. 



the stars should rotate close to the rates expected for 
single stars in the absence of a binary companion. Unfor- 
tunately, such an endeavor is difficult due to the reduced 
geometric alignment probability and limited availability 
of long-duration photometric monitoring of M-dwarfs at 
rapid cadences, where the majority of wide-field variabil- 
ity surveys are carried out in visible bandpasses, often 
yielding poor signal to noise ratios for M-dwarfs, which 
are faint at these wavelengths. 

The longest-period M-dwarf EB with a published orbit 
is T-Lyrl-17236 (Devor et al. 2008), which has an 8.4 day 
orbital period and components of 0.68 and 0.52 M0. The 
radii are indeed consistent with the theoretical models, 
but the observational uncertainties are still quite large 
in the discovery paper. This object is currently being 
observed by the NASA Kepler mission, so a more precise 
solution should be forthcoming. 

In light of the difficulty of identifying longer-period 
EBs, several authors have attempted to test the activity 
hypothesis using the existing EB sample to search for cor- 
relations between e.g. orbital period or activity measures 
and radius infiation (jLopez- Morales! 120071 : iKraus et al.l 
|2011| ). While there does appear to be evidence for such 
a correlation above 0.35 M0, below this mass the present 
sample is limited by small number statistics and the 
narrow range of orbital periods spanned. The Kepler 
mission appears to have ide ntified a large nu mber of 
very long-period EB systems (jPrsa et al.l 120111 : see also 
iCoughlin et al.l 120111 where more detailed solutions us- 
ing the Kepler photometry are presented for low-mass 
objects from this sample), but we are not aware of pub- 
lished dynamical masses and mo del- independent radii for 
these at the present time, and many of the M-dwarfs (es- 
pecially below 0.35 M0) are secondaries in systems with 
mass ratios much less than unity, where detecting the 
secondary spectrum to measure radial velocities may be 
challenging. 

An alternative method to determine M-dwarf radii is 
by making interferometric measurements of angular di- 
ameters for (usually isolated) stars. Combined with 
trigonometric parallaxes and minor assumptions about 
the stellar limb darkening, the physical radius of the 
star can be inferred in a mo del- independent fashion. 
The difficulty with this method comes in estimating 
the stellar mass. Dynamical measurements are usually 
not available, so one has to resort to inferring mass 
from other measurable stellar properties, necessarily with 
larger uncertainties. Typically, absolute magnitude in 
near-infrared passbands is used, which is found to corre- 
late well with stellar mass and has a relatively sma l l scat - 
ter. The empirical polynomials of iDelfosse et al.l (|2000[ ) 
are the most widely used, although there is some argu- 
ment in the literature as to how accurately these predict 
the mass, especially over a range of metallicity and spec- 
tral type (the uncertainties are probably > 10%). These 
concerns combined with small sample sizes mean inter- 
ferometry resul ts have thus far remained somewhat am- 
biguous (e.g. Dcm orv et all 120091: iBovaiian et al.l 1201 IL 
and references therein). 

In this contribution, we present the discovery of 
eclipses in LSPM J1112-h7626, a nearby double-lined M- 
dwarf system with an extremely long orbital period (com- 
pared to existing well-characterized M-dwarf EBs) of ap- 
proximately 41 days. At this large orbital separation the 



LSPM J1112+7626 



3 



system is not ex pected to h ave experienced significant 
tidal effects (e.g. IZahn|[l977l ). and is predicted to be nei- 
ther synchronized nor circularized, as verified by the non- 
circular orbit and apparent non-synchronized spin and 
orbit for one of the components as we shall describe. In 
the absence of strong tidal effects, it is expected that the 
angular momentum evolution of each component has not 
been perturbed significantly by the binary companion, 
and therefore the stars should possess magnetic prop- 
erties more representative of single M-dwarfs. In this 
contribution, we proceed to derive a high-quality spec- 
troscopic orbit and a photometric solution for the system. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Initial detection 

LSPM J1112-h7626 was targeted as part of 
routine operation of th e ME a rth transit survey 
(|Nutzman fc Charbonneaul l2008t llrwin et al.l l2009al) 
during the 2009-2010 season, with the first observations 
taken on UT 2009 December 2. Exposure times for 
the LSPM J1112-h7626 field were 27 s, taking three 
exposures per visit (the total exposure times are tailored 
for each target to achieve sensitivity to a particular 
planet size for the assumed stellar radius), with visits at 
a cadence of 20 minutes. 

Eclipses were first detected in LSPM J1112-H7626 on 
UT 2010 April 3 (this was a primary eclipse), and subse- 
quently confirmed on UT 2010 April 28th by observing 
the majority of a secondary eclipse egress. Both eclipses 
were detected by an experimental automated real-time 
eclipse detection and followup system we began to trial 
during the 2009-2010 season (this will be described in de- 
tail in a future publication when it has matured). Con- 
fident that the object was an eclipsing binary, but still 
with an ambiguous orbital period, we immediately com- 
menced radial velocity monitoring. The radial velocities 
ultimately determined which of the possible orbital pe- 
riods was correct, in conjunction with a further eclipse 
(which was also a secondary) observed on UT 2010 June 
8 with MEarth and the Clay telescope (see §2.2|) . 

2.2. Secondary eclipses 

By the time the orbital ephemeris was well-known, the 
observing season on LSPM J1112-I-7626 from Arizona 
was ending, and the weather had begun to deteriorate 
due to the arrival of the monsoon. We therefore took 
advantage of the high declination of our target and com- 
menced an observing program to target the secondary 
eclipses using the Harvard University 0.4m Landon T. 
Clay telescope, which is located on the roof of the Science 
Center on the Harvard campus near Harvard Square, 
Cambridge, Massachusetts. This instrument is normally 
used for undergraduate teaching. It is equipped with an 
e2v CCD47-10 1024 x 1024 pixel back-illuminated CCD 
packaged by Apogee Instruments (model AP47), which 
was operated binned 2 x 2 for a final pixel scale of « 1'.'59 
per pixel to reduce overheads. 

The location of the telescope is extremely light pol- 
luted, and combined with the red colors of our target, we 
decided to use the /-band filter to mitigate the effects of 
sky background and atmospheric extinct ion. The filter 
available was built to the prescription of 'Be sseUl (|199O0 
using Schott RG9 glass to define the bandpass, and shall 



hereafter be referred to as /bcssoU- While this filter was 
intended to approximate the Cousins / passband (here- 
after 7(7 ), the combination of the glass filter and CCD 
quantum efficiency yields a system response with signif- 
icant sensitivity at very red wavelengths not present in 
the original Cousins passband, and thus is a poor ap- 
proximation to the Cousins filter for very red stars such 
as M-dwarfs. We shall return to this issue in ^4.41 

Data were successfully gathered using the Clay tele- 
scope for the secondary eclipses on UT 2010 June 8, Oc- 
tober 9, November 19, and December 30th. The entire 
eclipse was visible for all but the June 8th event, although 
varying fractions of these were lost due to clouds. Expo- 
sure times were adjusted to keep the peak counts in the 
target below 30 000 to avoid saturation, and ranged from 
30 — 60 s. The typical cadence was 50 s and the FWHM 
of the stellar images was 2 — 3 binned pixels. No guiding 
was available, so the target was re-centered manually ap- 
proximately every 20 exposures to attempt to keep the 
drift below 5 pixels. 

Data were reduced and light curves produced using 
standard procedures for CCD data, following the method 
outlined in Irwin ct al. (2007). 2-D bias frames were sub- 
tracted as there was some bias structure evident, and we 
used the overscan region to correct for any bias level vari- 
ations. Dark frames of exposure matched to the target 
frames were subtracted to correct for the (significant) 
dark current and hot pixels at the — 20C operating tem- 
perature of the CCD, and the data were flat-fielded us- 
ing twilight fiat-fields in the usual way. No defringing or 
absolute calibration of the photometry were attempted. 
Fringing was clearly visible in the images, but at a low 
amplitude (3% peak-to-peak) and should not have a sig- 
nificant effect on the photometry due to the stabilized 
telescope pointing. 

We present the full light curve data-set used in this 
paper in Table [U and times of minimum light for the 
eclipses where a clear minimum was seen in Tabled) 

2.3. Primary eclipse 

Using the timing of the single visit within the pri- 
mary eclipse of UT 2010 April 3 from MEarth (see gjjl 
combined with modeling of the secondary eclipses ( H2.2 1 
and radial velocities it was possible to predict the 
primary eclipse times with reasonable accuracy. Dur- 
ing 2010-2011 full eclipses were observable at reasonable 
airmass (< 2.0) from eastern Europe, Scandinavia, and 
western Asia, with the high declination of the target 
favoring northerly latitudes during winter to maximize 
dark hours at reasonable airmass. 

The primary eclipse was observed using the 40 cm tele- 
scope at Hankasalmi Observatory, Finland on UT 2011 
January 15, starting approximately 10 minutes before 
first contact, and continuing until well after last con- 
tact. Data were taken using a similar /Besseii filter to the 
secondary eclipse observations (see ^2.2p with a front- 
illuminated Kodak KAF-IOOIE CCD packaged by Santa 
Barbara Instrument Group (SBIG; model STL-IOOIE). 
This system response should be very similar to the one 
used for the secondary eclipses, but nonetheless, any 
bandpass mismatch is a potential concern for the accu- 
racy of the photometric solution when combining data 
taken with different instruments, and will be discussed 
in more detail in ^A.4l 
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TABLE 1 
Light curve data. 



ID=' 


Orbit*^ 


HJD.UTC 


/ 




A mag'^ 


FWHM^ 


Airmass 


MeridFlag'= 


X 


y 


CM' 






(days) 


(mag) 


(mag) 


(mag) 


(pix) 






(pix) 


(pix) 


(mag) 


C 


1 


2455355.573517 


-0.0224 


0.0073 


0.000 


2.10 


1.2814 





100.427 


378.595 


0.0000 


C 


1 


2455355.574050 


-0.0168 


0.0077 


-0.002 


2.19 


1.2820 





100.715 


378.567 


0.0000 


C 


1 


2455355.574547 


-0.0161 


0.0079 


-0.000 


2.18 


1.2827 





100.936 


378.442 


0.0000 


C 


1 


2455355.575091 


-0.0176 


0.0079 


-0.033 


2.26 


1.2833 





100.904 


378.255 


0.0000 


C 


1 


2455355.575624 


-0.0165 


0.0081 


-0.010 


2.20 


1.2840 





100.900 


378.023 


0.0000 



Note. — Table [T] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 

Data-set identifier. C: Clay telescope secondary eclipse (magnitudes for each eclipse / night were separately nor- 
malized to zero median) H: Hankasalmi primary eclipse (normalized to zero median as above) O: MEarth out of 
eclipse (magnitudes in the instrumental system) P: MEarth primary eclipse (715nm long-pass, normalized as above) 
S: MEarth secondary eclipse (normalized as above) 

^ Orbit number (from zero at Tq). Given only for the data-sets containing eclipses. These are the numbers used in 
the figures showing the light curves. 

Photometric zero-point correction applied by the differential photometry procedure. Note: this has already been 
included in the / column and is for reference purposes only (e.g. detecting frames with large light losses due to clouds). 

FWHM of the stellar images on the frame. Zero in cases where the FWHM could not be reliably estimated. 
° Flag for detecting "meridian flip" on German equatorial mountings. Used only for the MEarth data, where the flag 
is when the telescope was observing in the usual orientation for the east side and 1 for the west side of the meridian. 
Note that there is not a one-to-one relation with the sign of the hour angle because the MEarth mounts can track 
approximately 5 degrees across the meridian before needing to change sides of the pier, so this flag can sometimes be 
even for (small) positive hour angle. 

f For the MEarth data only. Gives the "common mode" interpolated to the Julian date of the exposure. Derived from 
the average differential magnitude of all the M-dwarfs observed by all 8 MEarth telescopes in a given time interval. 
This should be scaled and subtracted from the / column to correct for the suspected variations in the MEarth bandpass 
with humidity and temperature (see text). 



TABLE 2 
Times of minimum light. 



HJD.UTC 

(days) 


Uncertainty 


Cycle 


Eclipse 


Instrument 


2455577.272644 


7 


7 


Primary 


Hankasalmi 


2455355.627938 


41 


1 


Secondary 


Clay 


2455478.725045 


42 


4 


Secondary 


Clay 


2455519.757412 


31 


5 


Secondary 


Clay 


2455560.789627 


11 


6 


Secondary 


Clay 


2455519.757351 


16 


5 


Secondary 


MEarth 


2455683.887000 


20 


9 


Secondary 


MEarth 



Note. — Est i mate d using the method of 

IKwee fc van WoerdenI (|1956[ ). over an interval of ±10"^ 
in normalized orbital phase. These times of minimum are 
not used in our analysis (we flt for the ephemeris directly 
from the data) and are reported only for completeness. 

Exposure times of 60 s were used, yielding a cadence of 
approximately 2 minutes, and the target was re-centered 
every 10 exposures to keep the drift below ~ 50 pixels. 
Data were reduced using 2-D bias frames, and master 
dark and flat-fleld frames. Light curves were then derived 
using the same methods and software as for the other 
data-sets. 

2.4. Out of eclipse and additional secondary eclipses 

Following the discovery of atmospheric water vapor 
indu ced systematic eff ects in the MEarth photometry 
(see llrwin et all 1201 ID , a customized /-band filter was 
procured for the start of the 2010-2011 season, with a 
bandpass designed to eliminate the strong telluric wa- 



ter vapor absorption at wavelengths > 900 nm, by using 
an 895 nm (50% transmittance point) interference cutoff 
on the same RG715 glass substrate as the original filter 
(where the long-pass glass filter defines the blue end of 
the bandpass). This was predicted to reduce the effect 
to largely insignificant levels, and the resulting system 
response approxi mates the Cousins / bandpass (as tab- 
ulated by BesseIil lT990l : see later in this section). 

Observations of LSPM J1112-h7626 using this new fil- 
ter commenced on UT 2010 October 29, and continued 
untfl UT 2011 May 18. Exposures of 2 x 97 s were used at 
20 minute cadence for times out of eclipse. 5283 useful 
data points were obtained at airmasses < 2, after dis- 
carding 400 bad frames (pointing errors, severe light loss 
due to clouds, and large scatter in the zeropoint solution 
for the differential photometry). 

Partial secondary eclipses were obtained with the same 
MEarth telescope on UT 2010 November 19, UT 2011 
February 9, and UT 2011 May 2. Due to known system- 
atic effects caused by persistence in the MEarth CCDs, 
lead ing to offsets in m agnitude depending on the cadence 
(see lBerta et al.|[201lD . and to remove any residual from 
the out-of-eclipse variations, we allow these high-cadence 
observations to have different photometric zero-points 
than the out of eclipse data (and each other) in the final 
light curve analysis. 

The out-of-eclipse data (excluding the high-cadence 
eclipse observations) were searched for periodi c varia- 
tions using the method outlined in llrwin et al.l ([2011). 
Correction for the "common mode" effect was still per- 
formed, as we discovered during the season that replacing 
the filters had in fact increased the size (and changed the 
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Fig. 1.— "F-test periodogram" for the MEarth 2010-2011 sea- 
son out of eclipse observations, plotting the square root of the F- 
test statistic comparing the null hypothesis of no variation to the 
alternate hypothesis of a sinusoidal modulation, as a function of 
frequency (top) and the corresponding window function (bottom). 
The vertical line indicates the best- fitting period of 64.8 days. 

sign) of the systematic effect rather than ehminating it. 
This problem is stih under investigation, but we suspect 
it is due to humidity and/or temperature dependence of 
the filter bandpass, an effect also seen in the original 
Sloan Digital Sky Survey (SDSS) Photometric Telescope 
filters (Doi et alJ 2Qf(l). 

After correcting for the "common mode" effect, the 
MEarth data show a strong near-sinusoidal modulation 
with a 65 day period (see Figure [U the data are shown 
in Q Slid a peak-to-peak amplitude of approximately 
0.02 mag. Marginal evidence for a similar periodicity is 
found in the 2009-2010 data, although the phase is not 
consistent and the amplitude appears to be much smaller. 
It is unclear if the amplitude change results from the 
differing filter systems, but the change in phase (if real) 
may be indicative of evolution of the photospheric spot 
patterns giving rise to the modulations. 

An improvement resulting from the new MEarth fil- 
ter is a greatly reduced color term when standardizing 
the photometry onto the Cousins system. The scatter in 
this relation also appears to be significantly lower, which 
may also result in part from improvements made to our 
fiat-fielding strategy during the 2010-2011 season. The 
MEarth observing software automatically obtains mea- 
surements of equat orial photomet ric standard star fields 
from the catalog of lLandoItl ()1992D every night, and these 
are used to derive photometric zero-points and calibrated 
photometry. The typical scatter in these zero-point so- 
lutions on a photometric night is 0.02 mag. We find a 
typical atmospheric extinction of 0.05 mag/airmass from 
MEarth data, which has been assumed henceforth. 

A significant issue with the equatorial standard fields 
is they contain very few red stars suitable for defining 
the red end of the photometric response, and so color 
equations derived only from these equatorial fields may 



TABLE 3 

Summary of the photometric and astrometric 
properties of the lspm ,11112 + 7626 system. 



Parameter 


Value 


02000''''' 


ll'il2'"42^.35 


fcooo""'' 


-l-76°26'56'.'3 


fj,a COS 


0'.'131 yr-l 




-O'.'IOS yr-l 


9" 


15.618 ±0.010 




14.195 ±0.012 




12.027 ±0.012 


Ic" 


12.14 ±0.05 




10.591 ±0.023 


-ff2MASS'^ 


9.989 ±0.021 


■f^2MASs'^ 


9.728 ±0.017 



- Equi nox J2000.0, epoch 200 0.0. 

b From lLepine fc Sh ara (2005). 

= From^ SpSS Data Release 7. Note that these 

uncertainties do not account for variability, and 

are ther efore underestimates. 

d From dl 

We quote the combined uncertainties from the 
2MASS catalog, noting that the intrinsic variabil- 
ity of our target means in practice that these are 
underestimates . 

be significantly in error when attempting to standard- 
ize M-dwarf photometry. Thankfully, MEarth serendip- 
itously targeted a number of M-dwarfs with photome - 
try available on the Cousins system from iBessell ()1990D . 
We selected all such objects in the 2010-2011 MEarth 
database with data available on photometric nights for 
calibration. Two objects with clearly unreliable MEarth 
data were discarded, leaving 8 objects. The following 
color equation was derived from these stars, which span 
Ic - -?f2MASS colors of 2.2 - 4.2: 

(0.085 ±0.019) 

+ (0.0204 ± 0.0064) (/MEarth - i^2MASs) 

In the case of LSPM J1112-I-7626, we estimate Ic = 
12.14 ± 0.05 using this relation, where we have infiated 
the uncertainty to account for variability and any errors 
introduced by the varying MEarth bandpass as discussed 
above. The observed Ic' — J, Ic — H and Ic — K colors are 
later used to infer effective temperatures and luminosities 
for the members of the binary, and thus the distance to 
the system. 

2.5. Literature data: proper motion and photometry 

In Table [3l we summarize the photometric and astro- 
metric system properties gathered from literature data 
and our photometry. This star lies within the area cov- 
ered by the Sloan Extension for Galactic Understand- 
ing and Evolution (SEGUE) stripe 1260 photometry, so 
we retrieved the PSF magnitudes from Data Release 7 
(lAbazajian et al.ll2009f ). This object is fiagged as sat- 
urated in i-band (the image also shows a very unusual 
PSF and may be corrupted), and the u-band photometry 
of M-dwarfs in SDSS is corrupted by a well-documented 
red leak problem, so we omit these filters. 

The moderately high proper motion of our target al- 
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lows the contribution of any background objects to the 
system hght to be constrained using previous epochs of 
imaging. In Figure [2l we show imaging data with the 
position and size of the MEarth photometric aperture 
(which is representative of the aperture sizes used for ah 
the photometry) overlaid. The DSSl image indicates the 
aperture is unlikely to contain any significant flux from 
background objects, and the SDSS image (which has the 
highest angular resolution) does not appear elongated. 

2.6. Intermediate-resolution optical spectroscopy: 

spectral typing and metallicity constraints 

A single epoch intermediate-resolution {R ~ 1200) 
spectrum of our target was obtained using the FLWO 
1.5 m Tillinghast re flector and FAST spectrograph 
(jFabricant et al.lll998D on UT 2011 March 16. The de- 
fault setting of a 300 g/mm grating and 3" slit was used 
to obtain two exposures of 300 s. The spectra were re- 
duced u sing standa rd procedures for long slit spectra in 
IRAPEI (lTodvlll993n . using a HeNeAr arc exposure taken 
inbetween the two target exposures to ensure accurate 
wavelength calibration. The 2-D spectra were divided 
by a flat field before extraction to remove the majority 
of the fringing seen at the reddest wavelengths. Relative 
flux calibration was performed using the spectrophoto- 
metric standard BD+26 2606. We show the resulting 
spectra in Figure [31 

Usin g the Ti05 band index defined by iReid et all 
()1995[) . we find a composite spectral type of M3.7 (note 
that the uncertainty in this value is at least 0.5 subclass), 
and using The Hammer (jCovev et all 120071 ). we find 
M3-M4 (visual inspection indicates the observed spec- 
trimi is closer to M4). The spectral type and the mea- 
sured colors are in good agreement using the empirical 
colors of Galactic disk M-dwarfs from Lcggctt (1992) for 
the Johnson/Cousins/CIT system and iWest et al.l ()2005h 
for the SDSS/2MASS. 

It is also clear from Figure [3] that no Ha line is seen in 
the spectrum. Such behavior is typical of inactive field 
dwarfs at these spectral types when observed at low reso- 
lution, where the equivalent width of the Ha absorption 
line b ecomes close to zero at M4-M5 (e.g. iGizis et"al] 
I2OOI see their Figure 5, noting that the Ti05 index 
was 0.42 for LSPM J1112-h7626). The presence of Ha 
emission in M-dwarf spectra is usually indicative of high 
activity levels. 

The relative strengths of the TiO and CaH bands in 
optical spectr a have been used to identify and classify 
M-subdwarfs ^Gizislll997l ). and more generally proposed 
as a method for estimating metallicities of M-dwarfs (e.g. 
iWoolf fc WallersteinI [200611 . While detailed calibrations 
(e.g. using band indices) do not yet appear to be avail- 
able, we performed a visual compa r ison b etween Figure 
[3l and the spectra shown in iGizisI (|1997f ) . The strong 
TiO bands seen in our target indicate it is not extremely 
metal poor. We therefore estimate it has metallicity 
[Fe/H] > —1, and is probably closer to solar (or greater). 

2.7. High-resolution optical spectroscopy: radial 

velocities 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement writh 
the National Science Foundation. 
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Fig. 2. — Images of LSPM J1112-I-7626 centered on the position 
as measured from the MEarth data. The circle shows the approx- 
imate position and size of the 8" (radius) photometric aperture 
used to derive the MEarth light curves. Data are from the first 
and second epoch Palomar sky surveys as provided by the Dig- 
itized Sky Survey (DSS; top and center panels), SDSS, and the 
MEarth master image (bottom panel). The approximate epochs 
are shown on the images, and all four panels have the same center, 
scale and alignment on-sky. 

High-resolution spectroscopic observations were ob- 
tained using the TRES fiber-fed echelle spectrograph on 
the FLWO 1.5 m Tilhnghast reflector. We used the 
medium fiber (2'.'3 projected diameter), yielding a re- 
solving power of i? ~ 44 000. 

Observations commenced on UT 2010 May 1, and con- 
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Fig. 3. — Intermediate-resolution FAST spectra of 
LSPM J1112-I-7626 showing the TiO and CaH molecular 
bands typically used for M-dwarf spectral classification. Two 
exposures were obtained and are shown offset in the vertical 
direction for clarity. Overlaid in gray are the locations of the 
Ha line, and the CaH2, CaH3 and Ti05 band indices defined by 
IReid et al.l lHOOST) . 

tinued until UT 2011 April 13, with 43 epochs acquired. 
Total exposure times ranged from 45 to 80 minutes, bro- 
ken into sequences of 3 or 4 individual exposures. The 
spectra were extracted using a custom built pipeline de- 
signed to provide precise radial velocities for echelle spec- 
tr ographs. The procedu res are described in more detail 
in lBuchhave et all (|2010[ ). and will be summarized briefly 
below. 

In order to effectively remove cosmic rays, the sets of 
raw images were median combined. We used a flat field 
to trace the echelle orders and to correct the pixel to pixel 
variations in CCD response and fringing in the red-most 
orders, then extracted one-dimensi onal spectra us i ng the 
"optimal extraction" algorithm of iHewett et al.l (|19850 
(see also iHorni 1 1 9861 ). The scattered light in the two- 
dimensional raw image was determined and removed by 
masking out the signal in the echelle orders and fitting 
the inter-order light with a two-dimensional polynomial. 

Thorium-Argon (ThAr) calibration images were ob- 
tained through the science fiber before and after each 
stellar observation. The two calibration images were 
combined to form the basis for the fiducial wavelength 
calibration. TRES is not a vacuum spectrograph and is 
only temperature controlled to 0.1°C. Consequently, the 
radial velocity errors are dominated by shifts due to pres- 
sure, humidity and temperature variations. In order to 
successfully remove these large variations (> 1.5 km s~^), 
it is critical that the ThAr light travels along the same 
optical path as the stellar light and thus acts as an effec- 
tive proxy to remove these variations. We have therefore 
chosen to sandwich the stellar exposure with two ThAr 
frames instead of using the simultaneous ThAr fiber, 
which may not exactly describe the induced shifts in the 
science fiber and can also lead to bleeding of ThAr light 
into the science spectrum from the strong argon lines, es- 
pecially at redder wavelengths. The pairs of ThAr expo- 
sures were co-added to improve the signal to noise ratio, 
and centers of the ThAr lines found by fitting a Gaussian 
function to the line profiles. A two-dimensional fifth or- 
der Legendre polynomial was used to describe the wave- 



length solution. 

3. SPECTROSCOPIC ANALYSIS 

Radial velocities were obtained using the two- 
dimensional cross-correlation algorithm TODCOR 
(jZucker & Mazeh 1994), which uses templates matched 
to each component of a spectroscopic binary to simul- 
taneously derive the velocities of both stars, and an 
estimate of their light ratio in the spectral bandpass. 

We used a singl e epoch observation of Barnard's 
star (jBarnardl 119161 : also known as Gl 699) taken on 
UT 2008 October 20 as the template for the TODCOR 
analysis. Barnard's st ar has a spectral type of M4 
(jKirkpatrick et al.1 119911 ). which should be a reasonable 
match to both components. Correlations were performed 
using a wavelength range of 7063 to 720lA in order 41 
of the spectrum to derive the velocities, since this region 
contains strong molecular features (mostly from TiO), 
which are rich in radial velocity information. We as- 
sume a barycentric radial velocity of — 110.3±0.5 km s~^ 
for Barnard's star (where the stated uncertainty refiects 
our estimate of the systematic errors), derived from 
presently unpublished CfA Digital Speedometer mea- 
surements spanning 17 years. 

To derive the light ratio, we ran TODCOR using a range 
of fixed values of this quantity, seeking the maximum cor- 
relation over all epochs. In this way, consistency between 
the different epochs is enforced, in accordance with the 
expectation that the true value of this quantity should 
be (approximately) constant. This procedure improves 
the stability of the solution and minimizes systematic er- 
rors in the velocities. We find a best-fitting light ratio 
of L2/L1 = 0.528 ± 0.05, where we adopt a conservative 
estimate of the uncertainty. 

The radial velocities derived from our analysis are re- 
ported in Table H) We omit four velocities with peak 
correlation values below Cpcak = 0.5 from the table and 
the remainder of the analysis as the low correlation val- 
ues indicate these velocities are likely to be unreliable or 
to exhibit large scatter. The uncertainty in the radial 
velocities (estimated from our fitting procedure; see Sj4]) 
is 0.26 km s~^ for the primary, and 0.58 km s~^ for the 
secondary, with an additional systematic uncertainty, es- 
timated to be 0.5 km s"^ (see above), in the velocity 
zero-point system. 

We find no evidence for "peak pulling" effects close 
to the systemic velocity, which might indicate the need 
for additional rotational broadening of the template to 
matc h the target. Barnard 's star is a very slow rotator 
(e.g. [Benedict et al.l[T998[ ). so this and the high corre- 
lation values obtained without need for additional rota- 
tional broadening indicate both components of our target 
have low vsini, as expected. 

4. MODEL 

For EB systems showing moderate or high eccentrici- 
ties, the commonly exploited property of being able to 
separate the photometric and spectroscopic solutions no 
longer holds as the eccentricity affects both. Typically, 
the orbital period P, epoch of primary eclipse To, and 
e cos u! (where e is eccentricity and w is the argument of 
periastron) are better determined by the photometry (by 
eclipse timing; the latter is related to the departure of the 
secondary eclipses in phase from being exactly halfway 
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TABLE 4 
Radial velocity data. 



HJD.UTC^ 


Phase'' 




V2'^ 


C 

pea 


fdavs) 


(k 


III S~ ) 


(km s^"*") 




2455317.6481 


0.67269 


-49.45 


-78.49 


0.560 


2455320.6671 


0.74626 


-41.11 


-90.08 


0.635 


2455321.6714 


0.77074 


-38.58 


-93.77 


0.643 


2455341.6568 


1.25780 


-79.89 


-33.62 


0.586 


2455342.6784 


1.28270 


-79.29 


-34.71 


0.506 


2455343.6545 


1.30649 


-78.44 


-36.09 


0.541 


2455344.6632 


1.33107 


-78.21 


-36.76 


0.577 


2455345.6848 


1.35597 


-76.63 


-38.66 


0.642 


2455346.6622 


1.37979 


-75.23 


-40.34 


0.536 


2455347.6488 


1.40383 


-74.26 


-42.12 


0.654 


2455348.6433 


1.42807 


-72.38 


-44.19 


0.602 


2455366.6561 


1.86706 


-35.28 


-97.31 


0.627 


2455367.6852 


1.89214 


-37.23 


-96.02 


0.588 


2455368.6748 


1.91626 


-40.62 


-90.98 


0.670 


2455369.6903 


1.94101 


-44.54 


-84.52 


0.576 


2455373.6462 


2.03742 


-65.72 


-53.84 


0.547 


2455374.6460 


2.06178 


-69.12 


-48.82 


0.662 


2455375.6455 


2.08614 


-72.87 


-43.53 


0.582 


2455581.9638 


7.11433 


-75.90 


-39.72 


0.750 


2455582.9270 


7.13780 


-77.72 


-37.42 


0.788 


2455583.8762 


7.16093 


-79.01 


-35.40 


0.811 


2455584.9160 


7.18627 


-79.99 


-34.08 


0.806 


2455586.9587 


7.23606 


-80.51 


-33.75 


0.854 


2455587.9329 


7.25980 


-80.20 


-33.98 


0.829 


2455605.9077 


7.69786 


-46.45 


-82.16 


0.816 


2455606.9374 


7.72296 


-43.79 


-87.07 


0.624 


2455607.9274 


7.74709 


-41.08 


-90.00 


0.797 


2455608.8926 


7.77061 


-38.34 


-93.22 


0.757 


2455611.8887 


7.84363 


-34.67 


-99.11 


0.731 


2455614.8719 


7.91633 


-40.32 


-91.05 


0.741 


2455615.8710 


7.94068 


-44.85 


-84.66 


0.822 


2455616.9096 


7.96599 


-50.16 


-76.85 


0.742 


2455617.8577 


7.98910 


-55.18 


-69.76 


0.788 


2455650.7940 


8.79179 


-36.93 


-95.63 


0.780 


2455651.7573 


8.81526 


-35.56 


-97.88 


0.732 


2455652.8008 


8.84070 


-34.89 


-98.32 


0.808 


2455656.7773 


8.93761 


-44.06 


-85.15 


0.814 


2455659.7417 


9.00985 


-59.54 


-62.78 


0.815 


2455664.7613 


9.13218 


-77.22 


-38.03 


0.826 



^ Heliocentric Julian Date of mid-exposure, in the 

UTC time-system. 

^ Normalized orbital phase. 

Barycentric radial velocity of stars 1 and 2. 
Peak normalized cross-correlation from TOD- 
COR. These are later used for weighting the radial 
velocity points in the solution. 

between the primary eclipses, although note that the 
commonly-used one-to-one relation between these quan- 
tities is approximate), and the esinw component is usu- 
ally best constrained by the velocities (but also affects 
the relative durations of the two eclipses, and in gen- 
eral eccentricity also influences their shape). Therefore, 
in order to leverage the best possible constraints on the 
system orbit, we adopt a joint solution of all the light 
curves and the radial velocities simultaneously. This also 
greatly simplifies the error analysis. 

We proceed by first discussing the basic model and as- 
sumptions in t ^4.1[ the method of solution in m.2l and 
then the effects of four potentially important sources 
of systematic error in the system parameters we derive: 
spots f !j4.3l) . bandpass mismatch f ji4.4p . limb darkening 
(t ^4.5|) . and third light (t ^4.6|) . In m.ll we summarize our 
best estimates of the system parameters and their corre- 
sponding uncertainties. 



4.1. Basic model and assumptions 

To model the system, we used the l ight curve genera- 
tor fro m the popular jktebop code (jSouthworth et aP 
l2004allb[ ). which is in turn based on the eclips ing bi- 
nary orbit program fEBOP: iPopper fc EtzellHOSl . This 
model approximates t he stars using biaxial ellipsoids 
(jNelson fc Davis! Il972t ) , which is appropriate for well- 
detached systems such as the present one. We modified 
the model to perform the computations in double preci- 
sion, to account for light travel time, which is necessary 
due to the large semimajor axis, and to incorporate a 
simple prescri ption for the effect of photospheric spots 
(discussed in !j4.3l) . The light curve model was coupled 
with a standard Keplerian orbital model for the radial ve- 
locities (also including the light travel time). We correct 
the orbital perio d from the heliocen tric frame to the sys- 
tem barycenter (|GrifRn et al.lll976l) when computing the 
physical parameters, following the usual convention of 
reporting the orbital period P as measured in the helio- 
centric frame (i.e. without the correction) in the tables. 

We fit aU of the Clay, Hankasalmi, and 2010-2011 sea- 
son MEarth data, the radial velocities, and the spectro- 
scopic light ratio simultaneously. We also include the 
single primary eclipse of UT 2010 April 3 to improve 
the constraint on the ephemeris, assuming the bandpass 
was approximately the same as the /bosscII filter (this 
makes little difference in practice due to the sparse cov- 
erage). The remainder of the MEarth 2009-2010 season 
data were omitted as these are of quite poor quality and 
were taken in a different bandpass from the remainder, 
so do not usefully constrain the model after accounting 
for the extra free parameters which must be added in 
order to fit them. 

Table [5] summarizes the assumptions made in our light 
curve model. The light curve parameter set generally fol- 
lows JKTEBOP, and has been chosen to minimize correla- 
tions between parameters where possible. We add a term 
ki{X — 1) to the model magnitudes, where X is airmass, 
and I is an index for the different light curves, to ac- 
count for differential atmospheric extinction between the 
target and comparison stars, presumably due to color- 
dependence. A different ki value was allowed for each 
light curve I. This has only been used for fitting the Clay 
secondary eclipse curves and 2010-2011 MEarth data (see 
i 32.2|) . where it seems to be necessary to include this term 
to obtain an acceptable fit to the out of eclipse portions. 
We are unsure as to the exact origin of the baseline cur- 
vature in the Clay data, and note that the choice of a 
linear function of airmass is arbitrary and may not be 
correct. This is further evidenced by the inconsistent 
sign and magnitude of this term in Table [6l We suggest 
the best approach to rectifying this problem would be to 
obtain higher quality secondary eclipse observations. 

Gravity darkening coefficients were fi xed at value s ap- 
propriate for convective atmospheres (jLucvl Il967f) , al- 
though gravity darkening is unimportant with the preci- 
sion of the present light curves at such low rotation rates. 
The reflection effect was computed rather than fitting for 
it, although again the large semimajor axis means it is 
unimportant. 

We used a square root limb darkening law 
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TABLE 5 

Parameters and priors used in the basic model. 



Parameter 


Value 


Prior^ 


Description 


i^MEarth 


varied 


uniform 


Central surface brightness ratio (secondary/primary) in /MEarth 


^^Bcsscll 


varied 


uniform 


Central surface brightness ratio (secondary/primary) in /bosscII 


+i?2)/a 


varied 


uniform 


Total radius divided by semimajor axis 


R2/R1 


varied 


uniform 


Radius ratio 


cos i 


varied 


uniform (isotropic) 


Cosine of orbital inclination 


e cos ui 


varied 


uniform in e,LU 


Eccentricity times cosine of argument of pcriastron 


esinw 


varied 


uniform in e,Lj 


Eccentricity times sine of argument of pcriastron 


L2/L1 


0.528 ±0.05 


Gaussian 


-fMEarth-band light ratio 


ui 


-0.0765 


uniform 


Primary limb darkening linear coefficient 


u[ 


1.0369 


uniform 


Primary limb darkening sqrt coefficient'' 


U2 


0.0561 


uniform 


Secondary limb darkening linear coefficient'' 


«2 


0.9335 


uniform 


Secondary limb darkening sqrt coefficient'' 


/3i 


0.32 




Primary gravity darkening exponenf^ 


/92 


0.32 




Secondary gravity darkening exponenf^ 


Kl 


computed 




Primary reflection effect coefficient 


K2 


computed 




Secondary reflection effect coefficient 


L3 





uniform 


Third light divided by total light'' 









Tidal lead/lag angle 


5 


2° 




Integration ring size'' 



Fi 


varied" 


uniform 


Primary rotation parameter 


ai 


varied" 


uniform 


Primary out-of-eclipse sine coefficient 


bi 


varied" 


uniform 


Primary out-of-eclipse cosine coeflicient 


F2 


varied" 


uniform 


Secondary rotation parameter 


0.2 


varied" 


uniform 


Secondary out-of-eclipse sine coefficient 


b2 


varied" 


uniform 


Secondary out-of-eclipse cosine coefficient 


q 


varied 


uniform 


Mass ratio 


K1 + K2 


varied 


modified Jeffreys 


Sum of radial velocity semiamplitudes 






Ka = 0.02 km S-' 


Prior P{Ki + K2) oc l/{Ki + K2 + Ka) 


7 


varied 


uniform 


Systemic radial velocity 


p 


varied 


uniform 


Orbital period (heliocentric) 


To 


varied 


uniform 


Epoch of primary conjunction (HJD_UTC) 


^SCC 


computed 




Epoch of secondary conjunction (HJD_UTC) 




varied 


uniform 


Baseline magnitude for light curve 


si 


varied 


Jeffreys 


Error scaling factor for light curve 


kl 


varied 


uniform 


Airmass coeflicient for light curve I (where used)' 


Ci 


varied 


uniform 


Common mode coeflicient for light curve I (where used) 


Sr 


varied 


Jeffreys 


Error (if Cpcak = 1) for radial velocity curve r' 



See §4.21 for discussfon. 
'' Varied in later sections. 

" Values appropriate for convective atmospheres. 

'' Our tests indicate this value results in a maximum error of 3 x 10~^ mag in the eclipse depth 
(compared to a model with 0.1°), which should be sufficient for our photometry. 
" Parameters varied only for one of the stars in each simulation. 

' I and r are indices for the different light curves and radial velocities for each component of the 
system. 



(|Diaz-Cordoves fc Gimene3ll992[ ). of the form: 



1 - u(l - /i) - u{l - ^) 



(1) 



where /a(m) is the specific intensity, and ^ = cos6', where 
9 is the angle between t he surface normal and the line 
of sight. As discussed bv lvan Hammd (|1993[ ). the square 
root law is a better approximation to the specific inten- 
sity distribution given by model atmospheres of late-type 
stars in the NIR than the other common two-parameter 
laws (quadratic and logarithmic) implemented in jkte- 
BOP. We also verified this by comparing the iClaretl 
(|200Clf) four-parameter law (which we assumed to be the 



best representation of the original PHOENIX model) with 
the two-parameter laws for typical M-dwarf stellar pa- 
rameters. 

For our baseline model, we fixed the limb darkening co- 
efficients to values for the Ig filter d erived from phoenix 
model atmospheres bv lClaretl (| 200Cf). using Toff — 3130 K 
for the primary, and Teff = 3010 K for the secondary, 
log (7 = 5.0 and solar metallicity. We assume the same 
limb darkening law applies to all the / filters (the error 
introduced by this assumption is examined in §4.51) . 

No evidence was found for any background sources in 
the photometric aperture (see i 32.5p , and the radial veloc- 
ity observations do not show additional lines or velocity 
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drifts that might indicate a comoving, physically associ- 
ated companion to the EB system. We therefore assume 
no third light (L3 = 0), but see also 



4.2. Method of solution 

To derive parameters and error estimates, we adopt 
a variant of the popular Markov Chain Monte Carlo 
(MCMC) analysis frequently applied in cosmology and 
for analysis of exoplanet radial velocity and light curve 
data (e.g. |Ttegmark et al.' '2004'; 'Ford 20051). We use 
Adaptive Metropolis (Haario ct al. 2001) rather than the 
standard Metropolis-Hastings method with the Gibbs 
sampler. The adaptive method allows adjustment of 
the proposal distributions while the chain runs, and is 
thus simpler and faster to operate, obviating the need 
for manual tuning of the proposal distributions typical 
of the more conventional methods. This method uses a 
multivariate Gaussian proposal distribution to perturb 
all parameters simultaneously. 

Compared to the simpler Monte Carlo and bootstrap- 
ping methods (e.g. as implemented in jktebop) these 
methods have an important advantage for the present 
case where a heterogeneous set of light curves and ra- 
dial velocities must be analyzed: it is possible to esti- 
mate the appropriate inflation of the observational un- 
certainties self-consistently a nd simultaneo usly with the 
other model parameters (e.g. lGregorvl2005D . This allows 
the correct relative weighting of the various observational 
data-sets to be decided essentially by their residuals from 
the best-fitting model, and the uncertainty in this weight- 
ing to be propagated to the final parameters. Of course, 
this assumes the correct model has been chosen for the 
data, so these methods must be used carefully. 

In the Monte Carlo simulations, we used the standard 
Metropolis-Hastings acceptance criterion, accepting the 
new point with probability: 



p — 



l,P(M,+i|D)>P(M,|D) 
n^^±im,P(M,+,\D)<PmD) 



where P{Mi\D) and P{Mi+i\D) are the posterior prob- 
abilities of the previous and current points, respectively, 
Ali represents the zth model, and D the data. The previ- 
ous point in the chain was repeated if the new point was 
not accepted (this is required to correctly implement the 
method). The appropriate ratio of posterior probabilities 
("Bayes factor") was: 



P(M,+i|Z?) _ P(M,+i) 



P{M,\D) 



P{Mi 



exp(-AxV2) 



(2) 



where the factors in the right-hand side are the ratio of 
prior probabilities, and the usual (least squares) likeli- 
hood ratio, respectively. 

Priors assumed for each parameter are detailed in Ta- 
ble [SJ and are chosen to be uninformative. In most cases, 
this choice is a uniform improper prior (labeled "uni- 
form" in the table) . For the orbital inclination, we adopt 
a uniform prior on cos i, which produces an isotropic dis- 
tribution of orbit normals, and for the eccentricity, we 
use a uniform prior on e and uj, rather than on the jump 
parameters e cos w and e sin u) (the la tter would produce 
a prior in e,a; pro portiona l to e; e.g. lFordll2006l) . 

As discussed bv lGregorvl ()2005D . the appropriate choice 
of uninformative prior for "scale parameters" such as ra- 



dial velocity amplitude or orbital period is a Jeffreys prior 
(or a modified Jeffreys prior for parameters which can be 
zero). This enforces scale invariance (equal probability 
per decade). In the present case, some of these parame- 
ters (particularly the period) are so well-determined the 
choice of prior is unimportant, so for simplicity we have 
only used non-uniform priors on the less well-determined 
parameters where the prior i s important. 

As discussed earlier, and in lGregoivl ()2005[ ). it is possi- 
ble to fit for the appropriate inflation of the observational 
uncertainties simultaneously with the model parameters. 
This has been done via the Sk parameters shown in Table 
[HI It is conventional to apply this inflation by multipli- 
cation for light curves, and for RV by adding an extra 
"jitter" contribution in quadrature. We follow these con- 
ventions here. Note that additional factors appear in the 
likelihood ratio of Eq. ([2|) when doing this, which have 
been subsumed into the Ax^ in the equation. 

For the RV, it is difficult to estimate reliable uncer- 
tainties from cross correlation functions, so instead we 
derive them during the simulation. We weight each 
point by C^^^y^ (the square of the peak normalized cross- 
correlation value; see Table |4]) , where we find Cpeak is 
approximately proportional to the signal to noise ratio in 
typical observations at similar signal to noise and peak 
correlation values to those seen in this work. This proce- 
dure is essentially equivalent to photon- weighting, except 
the cross-correlation accounts for all sources of uncer- 
tainty rather than only photon noise. The uncertainty 
corresponding to a peak correlation of unity is derived 
during the simulation, and we allow separate values of 
this quantity for each star, labeled si for the primary, 
and S2 for the secondary, in the tables. The resulting 
per-data-point uncertainties used in the fitting are given 
by (Jij = Si/Cpcak.j for the ith star and j'th radial velocity 
point. 

All Monte Carlo simulation results reproduced in this 
paper are derived from chains of 2 x 10® points in length, 
where the first 10% of the points were used to initial- 
ize the parameter covariance matrix for the Adaptive 
Metropolis method (starting from the initial parame- 
ters and covariance matrix derived using a Levenberg- 
Marquardt minimizatiorQ), and the next 40% were used 
to "burn in" the chain. This was found to be sufficient 
to ensure the chains were very well converged, and all of 
these first 50% of the points were discarded, leaving the 
remainder (10® points) for parameter estimation. Corre- 
lation lengths in all parameters were < 200 points. 

We report the median and 68.3% central confidence 
intervals as the central value and uncertainty for all pa- 
rameters. Reduced values for all the model fits were 
unity, by construction. 

4.3. Spots 

In LSPM J1112+7626, the presence of spots is clearly 
indicated by the observed out of eclipse modulation de- 
scribed in W2.A\ 

Spots present severe difficulty for the analysis of eclips- 
ing binary light curves because they cause systematic 
errors in the measurements derived from the eclipses, 
principally their depths, which determine the geomet- 

* |http: //www. Ics .forth.gr/~lourakls/levmar/ 1 
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ric parameters J, + R2)/a and i (and by exten- 
sion, R2/R1). Spots occulted during eclipse temporar- 
ily reduce the depth as they are crossed, and spots not 
occulted during eclipse increase it, because the surface 
brightness of the photosphere under the eclipse chord is 
greater than would be inferred from the out-of-eclipse 
baseline level. 

The difficulty in solving for physical system param- 
eters arises because the true spot distribution is not 
usually known. Out-of-eclipse modulations, and in non- 
synchronized systems, modulation of the eclipses them- 
selves, are sensitive only to the longitudinal inhomogene- 
ity in the spot distribution, except in special cases where 
spots are crossed during eclipse and cause detectable de- 
viations from the usual light curve morphology (this can 
be difficult to detect in systems with radius ratios close to 
unity, because the deviations then last a large fraction of 
the eclipse duration). Any longitudinally homogeneous 
component, such as a polar spot viewed equator-on, or 
a homogeneous surface coverage of small spots, cannot 
usually be detected from light curves. 

This can cause undetected systematic errors in the 
radii an d effective temperat ures derived from the light 
curves. [Morales et al.l (|2010f ) discuss this issue in some 
detail for the best-measured literature EBs, finding this 
effect could produce up to 6% systematic errors in the 
radii for stars with spots of 30% filling factor, when the 
spots are concentrated at the pole. 

The influence of the spectroscopic light ratio is not 
commonly discussed with regard to spots, but this is im- 
portant because it provides complementary information 
to the photometry, on the difference in spot coverage 
between the two stars. For example, consider the case 
of a large coverage of longitudinally homogeneous, po- 
lar spots that are not crossed during eclipse. If these 
spots are distributed on both stars, the light ratio is un- 
affected but both eclipses will appear deeper. However, 
if the spots are on only one star, the effect on the pho- 
tometry is the same, but the light ratio is altered be- 
cause the spotted star appears darker. This will produce 
a discrepancy between the light ratio derived only from 
the light curve parameters, and the spectroscopic value. 
More generally, it is possible to place constraints on the 
difference in the overall spot coverage between the binary 
components, using the spectroscopic information. 

One unusual feature of the present system among low- 
mass eclipsing binaries showing signs of spots is the spin 
and orbit appear to not be synchronized. Because of this, 
it is possible to measure eclipses at different rotational 
phases, and thus different spottedness of the visible stel- 
lar hemispheres. While this does not necessarily resolve 
the difficulty of determining the longitudinally homoge- 
neous component of the spot distribution, it docs pro- 
vide additional information on the longitudinally inho- 
mogeneous component, specifically which star the spots 
are located on. This information is not usually avail- 
able if only out of eclipse modulations are seen. It is 
common to assume the spots are on both binary com- 
ponents, which may be reasonable for near equal mass 
systems at short periods where tidal synchronization is 
expected to have occurred, but this is not a reasonable 
assumption for LSPM J1112-I-7626, where only one out- 
of-eclipse modulation is see n. The o bserved rotational 
evolution of M-dwarfs (e.g. Ilrwin et al.,,2011.) indicates 



it is extremely unlikely the two binary components could 
have the same rotation period (and phase!) by chance in 
the absence of tidal effects, which are not expected. 

Unfortunately, at the time of writing, only a single 
primary eclipse is available, so we are unable to take 
full advantage of the non-synchronized spin and orbit 
at present. Also, while multiple secondary eclipses were 
observed, and residuals from our best-fitting model (as- 
suming no changes in eclipse shape or depth) are seen, 
it is not clear if many of these are simply the result of 
systematic errors in the photometry, given that equally 
large deviations are seen out of eclipse. Constraining the 
spot properties by this method will be an important area 
for future work, and observing in multiple bandpasses 
would be advantageous as it may provide information on 
the spot temperatures. 

4.3.1. Spot model 

Conventionally, a Roche model such as the one im- 
plemented in the popula r Wilson-Devinney program 
(jWilson fc Devinnevl I1971D would be used to model a 
system with spots as these usually include a circular 
spot model and perform the necessary surface integrals 
over the stars; however the treatment of proximity ef- 
fects is completely unnecessary in the present system, 
and these models are extremely computationally inten- 
sive, especially for systems with eccentric orbits, which 
would make a detailed Monte Carlo based error analysis 
prohibitive. 

It is also not clear that the circular spot model 
with a small number of spots is realistic for late- 
type dwarfs. Observed light curves very rarely show 
the characteristic "eclipse-like" features with fiat base- 
line at maximum light, as would be predicted for a 
small number of near-equatorial spots. This indi- 
cates either tha t these objects hav e very large, po- 
lar spots (e.g. iRodono et al.l I1986D . such that some 
of the spot is always in view to the observer to pro- 
duce the continuous photometric modulations seen in 
light curves, or that t he surfaces of these objects have 
many small spots ( e.g. IBarnes fc Collier Cameronll2001l : 
IBarnes et al.l I20M) with the photometric modulations 
arising from a longitudinal inhomogeneity of the spot 
distribution. Modeling an ensemble of spots, as in the 
latter case, in detail would be prohibitive, as even single 
spot models are usually degenerate given limited light 
curve data. In our case, the degeneracies in fitting spot 
models would be further exacerbated by the availability 
of only single-band light curve information, meaning spot 
temperature and size would be essentially degenerate. 

Given these difficulties, we take an alternative ap- 
proach for the solution presented in this paper. Rather 
than trying to model the unknown spot distribution in 
detail, we simply attempt to incorporate the effect of 
spots into the uncertainties on the final parameters. We 
consider two models: (a) spots on the non-eclipsed part 
of the photosphere on both stars; and (b), a case in- 
tended to approximate the effect of spots covering the 
whole photosphere, again for both stars. In both cases, 
the models are required to reproduce the observed out of 
eclipse modulation. It is necessary to run two models in 
each spotted case because it is not known which star the 
out of eclipse modulation originates from. 

We assume a simple sinusoidal form for the modula- 
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tions, which appears to be an adequate description of the 
available light curve data. The functional form adopted 
was: 



Gi sm 



f2'KF,t 



+6i cos 



/ 2TTF,t 



(3) 



where t is time, Li is the light from star i, and hi 
are constants expressing the amplitude and phase of the 
modulation, and Fi is the "rotation parameter" , the ratio 
of the rotation frequency to the orbital frequency. This 
expression is constructed to yield AL, = at maximum 
light, corresponding to the conventional assumption that 
no spots are visible when the star is brightest, which gives 
the minimum spot coverage necessary to reproduce the 
observed modulations. 

JKTEBOP computes the final system light (and thus the 
change in magnitude) by summing the out of eclipse light 
and then subtracting the eclipsed light. By modulating 
only the out of eclipse light in this summation using Eq. 
([2]), hypothesis (a) can be implemented, and hypothesis 
(b) can be implemented by also modulating the eclipsed 
light (changing the surface brightness under the eclipse 
chord) when the eclipsed star is star i. 

4.3.2. Results 

The effect of applying our spot model on the eclipses 
is shown in Figures |4] and [5l which were calculated us- 
ing the parameters of LSPM Jl 112+7626, comparing a 
model with spots to a model with identical physical pa- 
rameters, but without spots. This figure demonstrates 
that observing multiple eclipses at high precision would 
allow the star hosting the spots to be identified. 

We fit all four spotted models to the full data-set 
for LSPM J1112-I-7626 using the procedures already de- 
scribed. The results are given in Tabled and Figures [6l 
[7| and [8] show the data with a representative model, cor- 
responding to hypothesis (a), the "non-eclipsed spots" 
model, on the primary, overplotted. 

The 2010 November 19 partial secondary eclipse from 
MEarth (numbered 5 in the figure) provides the most 
discriminating power between the various models, par- 
ticularly taken in conjunction with the 2011 May 2 event 
(numbered 9), which dominates the fit. As shown in Fig- 
ure [71 these were taken at opposite extremes of the out 
of eclipse variation, with the 2010 November 19 event 
at maximum light and the 2011 May 2 event close to 
minimum light. Examining the sizes of the s parame- 
ter in Table El the "eclipsed spots" model on the sec- 
ondary, which is the only model presenting a significant 
detectable deviation in the secondary eclipses (as shown 
in Figure [5]) is disfavored by the data (see Figure |9]), hav- 
ing an SMEaith-20101119 valuc larger by almost Scr than 
the "non-eclipsed spots" model on the secondary. In- 
deed, the observations indicate that this eclipse needs to 
be made as shallow as possible in the model in order to 
produce a good fit, possibly even slightly shallower than 
the "non-eclipsed spots" model is able to produce (the 
discrepancy may arise from the imperfect modeling of the 
out-of-eclipse variation by a sinusoid) . The depth of this 
eclipse therefore argues that if the spots causing the out 
of eclipse modulation are on the secondary, they should 
be at latitudes not crossed by the eclipse chord. 

It is interesting to note that, because the radius ratio is 
quite close to unity, the range of allowed combinations of 




Fig. 4. — Predicted deviation in the primary (top) and secondary 
(bottom) eclipses resulting from spots on the primary star. The 
solid curves correspond to hypothesis (a) where the spots are not 
eclipsed, and the dashed curves to (b) where they are eclipsed. 
The two sets of curves are identical in the lower panel. The cycle 
numbers (integer part of the normalized orbital phase) are shown 
on the right of each panel. 

spot latitude and inclination angle for the secondary that 
would explain the observed out of eclipse modulations 
is therefore quite narrow, compared to the much wider 
range of possibilities on the primary that are still allowed. 
This mildly argues in favor of spots on the primary as 
being the more likely explanation, purely from geometry. 

Nevertheless, with the presently-available data, the 
other models do not show any clear disagreement with 
the observations, so we consider (from a purely observa- 
tional point of view) that all three are equally likely at 
present. 

4.3.3. Effect of additional spots 

The model we have considered so far uses the minimum 
spot coverage necessary to reproduce the observed out- 
of-eclipse modulations. As noted, e.g. by iBarnes et al.l 
(2011) and other authors, studies of cool stars sensitive 
to the total spot coverage typically find much higher fill- 
ing factors than photometric spot models. Therefore it 
is likely these stars have a significant longitudinally ho- 
mogeneous spot component. 

In this section, we consider the effect of adding ad- 
ditional spots in a longitudinally homogeneous fashion 
where they would not be detectable through modula- 
tions in the photometry. For simplicity, we have placed 
these spots on the same star and at the same latitudes 
as the longitudinally inhomogeneous component, noting 
that we are predominantly interested in the effect on the 
component radii, where it does not matter which of the 
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TABLE 6 

Derived parameters and uncertainties for the four spot configurations with no additional spots. 



Parameter'' 



Non-eclipscd spots 
Primary Secondary 



Primary 



Eclipsed spots 



Secondary 



•.^MEa 



(Rl +R2)/n 
R2/R1 
cos i 

P (days) 

To (HJD.UTC) 

Fi 

ai (mag) 
bi (mag) 
F2 

a2 (mag) 
b2 (mag) 

sooe'' 

SMEarth-20101119 

«MEarth-20110209 

■5MEarth-20110502 

■^MEarth-Pri mary 

^Hankasalmi 

■5CIay-20100608 

*CIay-20101009 

*CIay-20101119 

*CIay-20101230 

^MEarth 

CMEarth 

'sCIay-20100608 

fcciay-20101009 

fcciay-20101119 

fcciay-20101230 

e cos ij 
e sin til 
<? 

{K1+K2) (km s-i) 
7 (km s"-'-)'^ 

51 (km s^^) 

52 (km s~^) 

i n 

e 

u (°) 
a (Rq) 
L2/L1 

(HJD.UTC) 

Ml (Mq) 

M2 (Mq) 

{R1+R2) (R©) 

Rl (Rq) 
R2 (Rq) 



0.7999 ± 0.0064 
0.8255 ± 0.0066 
n ni ^^/is+0 000039 

O.Uit>54»_(, (,Q(,Q42 

0.787 ±0.017 

fO. 000103 
-0.000112 

41.032363 ± 0.000024 
2455290.04618 ± 0.00018 



0.004645^ 



0.8004 ± 0.0066 
0.8259 ± 0.0069 

U.U10001_Q (JQ(JQ44 

0.788 ±0.018 
004652+°-'""'^''^ 
41.032364 ±0.000024 
2455290.04618 ± 0.00018 



0.7996 ± 0.0064 
0.8252 ±0.0066 

U.U1000U_g QQQQ42 

0.788 ± 0.017 
004646+° °°°^°^ 
41.032365 ±0.000024 
2455290.04617 ± 0.00018 



0.8423 ± 0.0047 
0.8637 ± 0.0043 
m5661+° °°°°25 

0.888 ±0.020 

HO. 000051 
-0.000058 

41.032336 ± 0.000025 
2455290.04637 ± 0.00018 



0.005005^ 



0.6332 ± 0.0019 




0.6331 ±0.0019 




0.00602 ± 0.00067 




0.00599 ± 0.00067 




0.00629 ± 0.00062 




0.00632 ± 0.00062 






0.6332 ± 0.0019 




0.6327 ± 0.0020 




0.0125 ±0.0014 




0.0098 ±0.0012 




0.0131 ±0.0013 




0.0109 ±0.0011 


1.469 ±0.014 


1.467 ±0.014 


1.467 ±0.014 


1.467 ±0.014 


1.191 ±0.058 


1.192 ±0.057 


1.191 ± 0.058 


1.498 ±0.086 


1.548 ± 0.087 


1.548 ±0.087 


1.549 ±0.087 


1.577 ±0.089 


1.194 ±0.055 


1.194 ±0.056 


1.193 ±0.055 


1.260 ±0.064 


1.40 ±0.36 


1.39 ±0.35 


1.42 ±0.36 


1.39 ±0.35 


1.135 ±0.049 


1.136 ±0.049 


1.137 ±0.049 


1.143 ±0.050 


1.022 ±0.046 


1.023 ± 0.046 


1.022 ± 0.047 


1.034 ±0.046 


1.457 ± 0.048 


1.455 ± 0.048 


1.456 ±0.050 


1.498 ±0.050 


1.177 ±0.038 


1.179 ±0.038 


1.178 ±0.038 


1.210 ±0.040 


1.367 ± 0.034 


1.367 ±0.034 


1.367 ±0.034 


1.357 ±0.033 


-0.00659 ± 0.00037 


-0.00658 ± 0.00037 


-0.00659 ± 0.00038 


-0.00686 ± 0.00039 


0.9834 ± 0.0080 


0.9841 ± 0.0080 


0.9842 ± 0.0080 


0.9868 ± 0.0080 


-0.0076 ± 0.0031 


-0.0074 ± 0.0031 


-0.0075 ± 0.0031 


-0.0204 ± 0.0032 


0.0136 ± 0.0017 


0.0135 ±0.0017 


0.0135 ±0.0017 


0.0184 ±0.0017 


0.0042 ± 0.0014 


0.0041 ± 0.0014 


0.0041 ± 0.0014 


0.0031 ±0.0014 


-0.0228 ± 0.0014 


-0.0228 ±0.0014 


-0.0229 ±0.0014 


-0.0244 ± 0.0014 


0.152408 ±0.000033 


0.152406 ± 0.000034 


0.152409 ± 0.000034 


0.152327 ±0.000035 


0.1823 ± 0.0011 


0.1824 ±0.0012 


0.1823 ±0.0012 


0.1852 ±0.0012 


0.6958 ± 0.0020 


0.6958 ± 0.0020 


0.6958 ± 0.0020 


0.6958 ± 0.0022 


55.602 ± 0.087 


55.602 ± 0.087 


55.603 ± 0.087 


55.601 ±0.091 


-61.071 ±0.027 ±0.5 


-61.071 ±0.027 ±0.5 


-61.070 ±0.027 ±0.5 


-61.068 ± 0.028 ±0.. 


0.125 ±0.015 


0.125 ±0.015 


0.125 ± 0.015 


0.134 ±0.017 


0.280 ± 0.033 


0.279 ±0.032 


0.279 ± 0.032 


0.290 ± 0.034 



89.7339ing^| 
0.23763 ±0.00086 
50.10 ±0.18 
43.825 ± 0.069 
0.485 ±0.025 
2455314.59517 ±0.00015 

0.3951 ± 0.0022 
0.2749 ± 0.0011 

n fi81A+f'0020 
0.6814_Q QP20 
n qsi t;+0.0028 

U.0010_Q QQ3Q 

3+0.0044 
-0.0043 



0.2999"* 



an 700C+O.OO66 

0.23770 ± 0.00088 
50.12 ±0.19 
43.824 ± 0.070 
0.474 ± 0.025 
2455314.59516 ±0.00015 

0.3950 ± 0.0022 
0.2749 ±0.0011 

u.b»i&_Q 0021 

n QO1/I+0.0028 
U.dUUi_Q 0Q44 



SQ 7QQC+0.0063 

0.23762 ± 0.00087 
50.10 ±0.19 
43.825 ±0.070 
0.486 ± 0.025 
2455314.59516 ± 0.00015 

0.3951 ±0.0022 
0.2749 ±0.0011 

u.b»i&_o 0021 

n QO1/I+0.0028 
U.3814_p Q(,3^ 

U.dUUi_Q QQ43 



SQ 71 Q9+0.0033 
■^^-0.0029 

0.23977 ± 0.00090 
50.56 ±0.19 
43.800 ± 0.073 
0.636 ±0.031 
2455314.59536 ± 0.00015 

0.3944 ± 0.0023 
0.2744 ±0.0012 

n QfiQK+0.0035 
n Q9oc:+0.0043 

U.OZZO_Q QQ42 



^ Parameter names as defined in Table [5] and in the text. 
The acronym OOE is used as shorthand for "out of ecHpse" . 

Errors quoted for the 7 velocity are the random error (scatter), followed by our estimated systematic error in 
the velocity zero-point, from the assumed barycentric velocity of Barnard's star. 
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Fig. 5. — Predicted deviation in the primary (top) and sec- 
ondary (bottom) eclipses resulting from spots on the secondary 
star. Curves as Figure |4] The two sets of curves are identical in 
the upper panel. 

possible locations is used for the longitudinally inliomo- 
geneous spot component. 

As discussed earlier, the spectroscopic light ratio is sen- 
sitive to differences in the overall spot coverage between 
the two stars in some cases. We have explored this limit, 
and the effect of adding varying quantities of spots on 
the component radii, by subtracting an extra term Ci in 
Eq. ([3]) representing the fraction of stellar light removed 
by the longitudinally homogeneous spot component. 

For single late-type active dwarfs, it is typical to find 
spot filling factors of 20 — 50% from obser vations sensi- 
tive t o the entire surface spot coverage (e.g. lO'Neal et al.l 
|2004[) , although we note the available information for M- 
dwarfs is very limited. Photometric observations of M- 
dwar fs indicate spot ternperat ure contrast of up to 10% 
(e.g. iRockenfeller et al.l 12001 who find 4 - 8%). This 
yields a decrement of approximately 10 — 20% in the I- 
band stellar light, assuming parameters appropriate to 
the present system, with the lower end of the range pre- 
sumably more typical for less active stars. We there- 
fore consider values of q — 0.1 and 0.2. The latter cor- 
responds approxima tely to the spot levels discussed by 
[Morales et al.l ()201C0 for the existing short-period eclips- 
ing binaries. The results are given in Tables [7] and [H 

Firstly, it is clear from Table [7] that adding additional 
"eclipsed spots" has no effect on the radii or light ratio 
and merely alters the surface brightness ratio (this means 
it also changes the effective temperature ratio). This is 
also true when placing these spots on the secondary, but 
we do not show results for this model since it is unable to 
reproduce the observed eclipse light curves, as discussed 
in gXl 



Adding "non-eclipsed spots" affects the light ratio, and 
both component radii, slightly increasing the sum of the 
radii, increasing Ri, and decreasing i?2- The light ra- 
tio behaves as expected, becoming closer to unity as 
the primary is darkened by the addition of spots. Both 
C2 — 0.2 models with "non-eclipsed spots" (and to some 
extent, the C2 = 0.1 model) produce a light ratio that is 
marginally discrepant with the spectroscopic value (see 
f|3] or Table [5]) , but the other models agree reasonably 
well within the errors. 

Given the lack of clearly-detected perturbations in the 
best quality eclipse curves, it is reasonable to favor the 
"non-eclipsed spots" models over the "eclipsed spots" 
models. Assuming the Ci = 0.1 model is typical, and 
the Ci = and 0.2 models represent the likely range 
of values, we estimate the final system parameters and 
their uncertainties by combining the posterior probabil- 
ity distributions from all six "non-eclipsed spots" models 
in this section (all three values of Cj). These parameters 
and estimated uncertainties are reported in ^4.71 

4.4. Bandpass mismatch 

A concern when combining light curves taken with dif- 
ferent instruments, is the effect of any difference in the 
photometric bandpasses. In the present case, where the 
primary and secondary eclipses were (by necessity) mea- 
sured using different instruments, this predominantly af- 
fects the ratio of the eclipse depths, and thus the param- 
eter J, the ratio of central surface brightnesses. 

In order to determine the approximate size of the band- 
pass mismatch, we measured the difference in magnitude 
between LSPM J1112-I-7626 and a nearby, much bluer 
comparison star, at ll'*12"12M5 -f 76°27'33'.'8 (posi- 
tion from 2MASS; this is the bright star seen near the 
upper left corner of Figure [2|). The latter star has 
(J — i^)2MASS — 0.27 and is thus probably in the F 
spectral class. It is non-variable at the precision of the 
MEarth data (rms scatter 0.003 mag). 

Our measured magnitude differences were 
A/(Hankasalmi) = 0.36, A/(Clay) = 0.35, and 
A/(MEarth) — 0.55. The excellent agreement between 
the two /bcsscII filters (especially in light of the observed 
out of eclipse variations of our target) indicates they are 
most likely an extremely good match. 

Assuming a simple linear scaling, we estimate the 
change in J corresponding to the mismatch between the 
two /bcsscII filters is approximately 1/10 of that between 
the /bgssoU and /MEarth filters, or 5J ~ 0.003. This is 
unimportant compared to the other sources of uncer- 
tainty (see Table [6]). 

4.5. Limb darkening 
4.5.1. Differences between /bcssgII and Ic 

We first examine the effect of assuming the limb dark- 
ening law is the same in the two / filters. To do this, 
we compare the fitting results using our baseline model 
to usin g the coefficients for the z' SDSS passband from 
Clarej ( 20041 ). where the /bcsscII band should lie inbe- 
tween these two extremes. 

Using the observed colors of LSPM J1112-h7626, we 
estimate that approximately 1/3 the difference between 
the Ic and z' results is appropriate for the error intro- 
duced by our assumption of the Ic limb darkening law 
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Phase Phase Phase 

Fig. 6. — Eclipse light curves used in the model, with the best fit overlaid (the model is the first of the assumptions in Table (6] but the 
other model light curves are almost identical in appearance). In the top panels, the model is shown with red dots at the epochs of the 
observations, and in the other panels it is plotted as a continuous curve. The figure has been divided horizontally into one column for each 
combination of eclipse and passband, with the primary eclipses (all in /bcsscIi) plotted on the left, the secondary eclipses in /jviEarth ™ 
the center, and in /bgsscII on the right. Each column contains three panels showing the raw differential photometry (top), the same after 
correction to flatten the out of eclipse baseline (center; these corrections are the airmass term for the Clay and MEarth secondary eclipses, 
and the "common-mode" term for the MEarth secondary eclipses) and the residual in the bottom panel. We have offset the different 
eclipses vertically for clarity, and the cycle number (integer part of the normalized orbital phase) is shown on the right. The vertical bars 
indicate the approximate locations of the first and last contact points. 



TABLE 7 

Derived parameters and uncertainties showing the effect of additional spots on the primary. 







Non-eclipsed spots 






Eclipsed spots 




Parameter 


ci = 


ci = 0.1 


ci = 0.2 


ci = 


ci = 0.1 


ci = 0.2 



-^MEarth 
■^Bcsscll 

(Ri + R2)/a 
R2/R1 



0.7999 ± 0.0064 
0.8255 ± 0.0066 

n ni c:c:/io-l-0. 000039 

0.787 ±0.017 

n04645+°-°"°l''3 



0.8302 ± 0.0066 
0.8583 ± 0.0068 
ni 56'^2+° '""'°33 

U.U10DJZ_() QQQQgg 

0.766 ±0.016 

005219+° '""'°*^ 



0.8724 ± 0.0074 
0.9042 ± 0.0079 

n ni cyyc-l-0. 000037 
U.UlOf '0_o. 000039 

0.756 ± 0.017 

0.00&861_o^000072 



0.7996 ± 0.0064 
0.8252 ± 0.0066 
015550+0-000041 

U.U1000U_Q 000042 

0.788 ± 0.017 

004646+0-000105 



0.7196 ± 0.0060 
0.7427 ± 0.0062 
015551+0-000042 

U.UiOOOi_o 000044 

0.788 ±0.018 

004646+0-000109 

U.UU4D4D_o 000j^j^g 



0.6396 ±0.0051 
0.6602 ± 0.0054 
015550+0-000041 

U.U1000U_Q 000042 

0.787 ±0.017 

004644+0-000107 

U.UU4D*4_Q^000j_j_2 



89.7339l»:00«^ 
0.485 ± 0.025 



89.7010lHr7 
0.531 ±0.025 



0.611 ±0.032 



89.7338l«;*'6? 
0.486 ± 0.025 



HQ 7QQO-I-0.0067 

0.486 ± 0.026 



89.7339ing^l 
0.486 ±0.025 



(R1+R2) (R©) 
Ri (R0) 
R2 (R©) 



.-1-0.0020 
-0.0020 

--I-0.0028 
-0.0030 
2999+0-0044 



0.68143 
0.3815j 



n fiS/t7+0-0020 
O.684/_0 0020 
n oovo-l-0.0026 
O.JSl-S^o 0029 

2969+0-00-12 
u.zjDy_0 0030 



O fiQn4+0-0020 r, -,01 c+0.0020 

O.69O4_0^OO2i U.6S15_0,oo2i 



QQQ4+0.0030 
U.ciyj^_0 0033 

2970+0-0045 
U.zy (U_o 0043 



O QO1/I-I-0.0028 
O.3S14_0 003^ 

•5001+0-0045 
U.dUUi_o 0043 



0.6815 
0.3814; 



+0.0021 
0.0021 

+0.0029 
0.0032 



3001+0-0047 
U.JUUl_o 0044 



0.6815 
0.3814 



+0.0020 
0.0021 

+0.0028 
0.0031 



3000+0-0045 
U.JUUU_o 0043 
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TABLE 8 

Derived parameters and uncertainties showing the effect of 
additional spots on the secondary. 







Non-eclipsed spots 




Parameter 


C2 = 


C2 = 0.1 


C2 = 0.2 



-^MEarth 
"^Bcssell 

{Ri+R2)/a 
R2/R1 
cos i 



0.8004 ±0.0066 
0.8259 ± 0.0069 

U.U10001_Q gggQ44 

0.788 ±0.018 

fO. 00011 
-0.00012 



0.00465^ 



0.8157 ±0.0075 
0.8409 ± 0.0077 
n ni Kcofi+o. 000043 

0.777 ± 0.016 

f-o. 00011 
-0.00012 



0.00493^ 



0.8271 ±0.0081 
0.8516 ±0.0084 

015601+'' °°°°'"^ 
U.UiODUi_o. 000045 

0.757 ±0.014 
u.uuoici_(, 00012 



L2/L1 



SQ vQqK+o.ooee 

o».loou_o 0Qg2 

0.474 ±0.025 



oQ 7-177+0.0086 
S9.7177_o.oo64 
0.423 ±0.022 



3+0.0067 
-0.0066 

0.360 ±0.017 



89.7063^ 



{R1+R2) (R©) 
Ri (R©) 
R2 (R©) 



0.6815 



+0.0021 
0.0021 

+0.0028 



0.3814_Q 0032 
3001+° "°*^ 



0.6829 



+0.0021 
0.0021 

+0.0026 



0-3844_ooo29 

0.2985_nnn/n 



0.6834 



+0.0022 
0.0022 
+0.0021 



0.3892_o 0024 

2942+° °°''° 
U.Z3'iZ_o 0037 
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Phase 

Fig. 8. — Radial velocity curves used in the model, with the best fit overlaid (top) and residual (bottom). The primary velocities are 
shown with square symbols, and the secondary velocities with round symbols. The residuals for the two stars have been offset for clarity. 



TABLE 9 

Comparison of derived parameters and uncertainties 
for ic and z' limb darkening laws using the 

"non-eclipsed spots" model on THE PRIMARY. 




Parameter 




z' 


■^MEarth 


0.7999 ± 0.0064 


0.7829 ± 0.0085 


--^Bessell 


0.8255 ± 0.0066 


0.8073 ± 0.0091 


(Ri+R2)/a 


ni554S+''-°°°°33 
U.Uia04»_g Q(,Q(,42 


m 53S7+°-°°°°53 


R2/R1 


0.787 ±0.017 


0.773 ±0.022 


cos i 


0.00464l°°««}; 


00440+° °°°'^ 


in 


89.7339t°°°«^ 


89 7481+0 0096 


L2/L1 


0.485 ± 0.025 


0.462 ±0.031 


iRi+R2) (R©) 


0.6814_g (,p20 


6744+°-°°25 


Ri (R©) 


sS15+° °°2* 

U.0010_g 


n Qsn7+0 0033 


R2 (Rq) 


2999+° °°"'' 
u.^aa»_g QQ43 


n 9QQQ+0.0059 
u.zaoo_(3 QQgg 



0.594 0.596 0.598 



Fig. 9. — As Figure|6l except for the "eclipsed spots" model on 
the secondary, showing the deviations in the MEarth 2010 Novem- 
ber 19 event. 



in fitting the /bcsscII data. Comparing the resuhs for 
the two bands in Table ^ we find this source of error in 
the radii is negUgible (ahhough it does affect J and i at 
close to la). While the radius sum error is comparable to 
the individual estimated uncertainties in the table, the 
uncertainty in this parameter for our adopted solution 
combining six spot configurations is larger, and the limb 
darkening contribution is then less important. 

4.5.2. Errors in the atmosphere models 

We have assumed limb darkening coefficients fixed to 
the theoretical values from phoenix model atmospheres 
throughout the analysis. We note that the same atmo- 
sphere models have issues reproducing the observed spec- 
tra of M-dwarfs in the optical, which raises the possibil- 
ity of systematic errors in the limb darkening law. In the 
literature, this is conventionally addressed by analyzing 
photometry in multiple passbands, which is sensitive to 
differences in the limb darkening as a function of wave- 
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length. We do not have multi-band photometry, but note 

this could be an important area fo r future work. 

It has long been recognized (e.g. iNelson fc Davislfl972l 
and references therein) that changes in the radius of the 
eclipsed star mimic changes in the hmb darkening. In 
a system with near-equal stars such as the present, it 
is reasonable to expect any systematic error in the limb 
darkening from the models should affect both stars sim- 
ilarly, and thus to first order, the dominant effect will 
be on the sum of the radii (and the inclination, but any 
uncertainties here have negligible effect on the final pa- 
rameters) . 

We verified this using simulations where the square 
root limb darkening coefficients were varied, finding a 
strong correlation between (i?i + R2)/a and the integral 
of the specific intensity over the stellar disc, which is 
1 — it/3 — u'/5 for the square root law (it is also corre- 
lated to a lesser extent with J and cosi). This quantity 
is the normalization term in the eclipse depths in the 
photometric model, so it is not surprising that it should 
directly infiuence the quantities derived from the abso- 
lute eclipse depth. This is predominantly a concern for 
the interpretation of the sum of the radii: we find even 
quite large changes in the limb darkening law do not 
significantly alter the individual radii compared to their 
uncertainties in our adopted model. 

4.6. Third light 

We now examine our assumption of no third light. 
While the proper motion evidence (see Figure [2]) argues 
it is unlikely there are any background stars contribut- 
ing at a significant level in the photometric aperture, the 
relatively low proper motion of our target does not allow 
this possibility to be completely eliminated given the lim- 
ited angular resolution of the first epoch imaging data. 
Common proper motion companions are also permitted 
at intermediate separations below the angular resolution 
of the SDSS data, but still in wide enough orbits or with 
mass ratios much less than unity, where they would not 
be detected by the presence of additional lines in the 
spectra or radial velocity drift over the approximately 1 
year baseline available. 

Therefore, to examine any constraints on third light 
which can be placed directly from the light curves, and 
the effect on the derived parameters, we ran an addi- 
tional set of Monte Carlo simulations using the basic 
"non-eclipsed spots" model on the primary discussed in 
i i4.3.2l and Table IH allowing L3 to vary. For simplicity, a 
uniform prior was assumed. We find L3 < .029 at 95% 
confidence. As discussed bv INelson fc DavisI ( |1972) . the 
dominant parameter mimic is between third light and 
inclination, and our results confirm this, finding the un- 
certainties in cos i and the central surface brightness ratio 
parameters were slightly inflated. The other light curve 
parameters and radii derived from these chains are in- 
distinguishable within the uncertainties from the results 
where L3 was not varied. 

4.7. Summary and adopted system parameters 

Of the sources of uncertainty we have considered, spots 
dominate. Since the stars hosting the spots and the spot 
configuration are mostly unknown, we adopt the union of 
the six spotted solutions discussed in ij4.3.31 giving them 



TABLE 10 

Adopted physical parameters and uncertainties, 
combining six spot configurations. 



Parameter 




Value 






COS'' 







Ml (Mq) 

Ah (Mq) 

(fli + R2) (R©) 
Ri (R©) 
R2 (R©) 

log 31 
log 92 

Teff.l^' (K) 
Teff,2" (K) 
T'cff,2/T'off,l 

-^^boI,l (Lq) 
-^^boI,2 (L©) 

(m - M) 
d (pc) 

C/t^ (km/s) 
V (km/s) 
W (km/s) 



0.3946 ± 0.0023 
0.2745 ±0.0012 

0.3860l°:»g 

n 9Q7Q+0-0049 

4.861 ±0.012 
4.929 ±0.014 



3061 ± 162 
2952 ± 163 

0.0119 ±0.0025 
0.0061 ± 0.0014 

3.40 ±0.22 
48.1 ±5.0 

56.1 ±3.1 
-44.4 ± 1.7 
-10.8 ±3.2 



3191 ± 164 
3079 ± 166 
960+° '''** 

U.9DU_Q 

0.0141 ± 0.0028 
0.0073 ± 0.0016 

3.53 ±0.22 
51.0 ±5.2 

57.6 ±3.2 
-44.8 ± 1.8 
-9.1 ± 3.4 



^ Effective temperatures and bolometric cor- 
recti ons used. COS: Casag rande et al. (2008); 
L96: iLeggett et all (|1996( ). A 150 K system- 
atic uncertainty was assumed for the effective 
temperatures, and the color equations given in 
the 2MASS explanatory supplement were used 
to convert the observed JHK photometry to 
the CIT system when using the L96 tabula- 
tion. 

^ Adopting the definition of positive U to- 
ward the Galactic center. Calcula t ed usi ng the 
method of IJohnson fc Soderbloml (|1987[) . 

equal weight. The final parameter estimates and their 
uncertainties are reported in Table 1101 derived from the 
posterior samples produced by the Monte Carlo simula- 
tions, adopting the median and 68.3% central confidence 
intervals as the central value and uncertainty (see ^A.2\ . 

We compute UVW space motions for 
LSPM J1112-H7626 usin g the method of 
IJohnson fc Soderbloml ()1987l ). with the distance 
from the eclipsing binary solution, 7 velocity from Table 
[Bl and position and proper motions from Table |31 We 
find this object is in the old Galactic disk pop ulation, 
following the method and definitions in ILeg gett! (jT992). 

For very precise work, even matters as seemingly trivial 
as physical constants can be important, so we briefly dis- 
cuss our assumptions in this regard. We adopt the 2009 
lAU values of GMq and the astronomical unit (see the 
Astronomical Almanac 2011), and follow Cox (2000|) in 
adopting the value of the solar photosp heric radius from 
iBrown fc Christensen-DalsgaardI ()1998r FI. For the solar 
effective teniperatu re, we use the value of 5781 K from 
iBessell et al.l ()1998l ). which is based on measurements of 

^ A useful compilation of solar data may be found 
in "Basic Astronomical Data for the Sun", by Eric Ma- 
majek (University of Rochester, NY, US A) , available at 
http : //www .pas . rochester . edu/ ~ emama j e k/sun . txt | 
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the solar constant by [Duncan eFall (fT98l) . 

The final constant, the solar bolometric absolute 
magnitude, is a m ore thorny issue, a s discu s sed b y 
iBessell et all (|1998l ). and extensively by iTorresI (|2010( ). 
When using tables of bolometric corrections, it is vital 
to adopt the consistent value of this quantity in accor- 
dance with the table. We therefore use the appropriate 
values in each case in the following section. 

4.8. Effective temperatures, bolometric luminosities and 

distance 

In order to estimate effective temperatures (and thus, 
bolometric luminosities, and the distance), an external 
estimate of the effective temperature of one of the bi- 
nary components is required, in addition to bolomet- 
ric corrections. For M-dwarfs, there are large sys- 
temat ic uncertainties in t he effective temperatur e scale 
(e.g. 'Al lard et all 119971: iL^man fc Rieki |1998(), with 
disagreement at the few hundred degrees Kelvin level 
among different authors. Many of the early difficulties 
were due to the lack of model atmospheres (which are 
usually needed to integrate the full SED from the avail- 
able measurements to obtain Lboi), and significant im- 
provement on this front was made in the 1990s. The 
availability of a much larger sample of angular diame- 
ter measurements from interferometry should further im- 
prove the situation as this provides a much more direct 
method to estimate T^a for single stars, but there are 
relatively few temperature scales available at the present 
time using this information. 

We show results from two different inversions to illus- 
trate the typical range of parameters. Both scales we 
use were derived with model atmospheres rather than 
blackbodies, and in both cases we use the measured 
Ic — J2MASS, Ic — -ff2MASs and Ic — -f'^2MASS colors 
from Table [3] in conjunction with the component radii 
and /c-band light ratio from the adopted eclipsing bi- 
nary solution. The results from these different pairings 
of bandpasses were found to be consistent within the un- 
certainties, so we report the union of the results. 

The two sets of effective te mperature and bolometri c 
corrections adopted were from iCasagrande et al.l ()2008D . 
which is a recent determination using interferometric an- 
gular diameter mea surements to derive Teg, and from 
iLeggett et al.l (|1996l ). which is the basis for several recent 
works an d tabulations of th e effec ti ve tempe r ature scale, 
includi ng 'Luhman fc Riekc' ()1998[ ). iLuhmanI ()1999[ ). and 
iKraus fc Hillenbrand (2007). 

For the iLeggett et al.l (| 19961 ) scale, we fit polynomi- 
als to their tabulated data, omitting three objects these 
authors found to be metal poor fr om our fits: Gl 129 , 
LHS 343, and LHS 377. We follow ILeggett et al.l dUp) 
in adopting a systematic uncertainty of 150 K, which 
appears to be consistent with the differences we find be- 
tween the two determinations. 

Our resuhs for LSPM J1112-I-7626 are reported in Ta- 
ble [lOl 

5. DISCUSSION 

We now compare our resul ts for LSPM J 1112- I-7626 
with theoretical models from iBaraffe et al.l (|1998[ ) and 
other literature objects with shorter periods summarized 
in Table 1111 We first show the conventional mass-radius 
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Fig. 10. — Mass-radius relation for components of detached, 
double-lined, double-eclipsing binary systems below 0.4 Mq. 
LSPM J1112-I-7626 is shown as bold, black points, and the lit- 
erature systems from Table 1111 are shown in gray, excluding 
CCDM J04404-(-3127B,C, which does not have a full solution avail- 
able at present, but including KOI 126. The lines show stellar evo- 
lution models from Baraffe et al. ( 1998) for 10 Gyr age, [M/H] = 
(solid line), [M/H] = -0.5 (dashed line), and the [M/H] = model 
with the radius inflated by 5%, corresponding to p = 1.05 (dotted 
line) . 
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Fig. 11. — Mass-effective temperature relation for components 
of detached, double-lined, double-eclipsi ng b inary systems below 
0.4 M0. Lines and points as Figure |lO] KOI 126B,C and 
IRXS J154727.5+450803 do not have reported effective temper- 
atures, so they do not appear on this diagram. The dotted line 
shows the prediction from assuming p = 1.05 and that the bolo- 
metric luminosity is preserved. 

diagram in Figure llOi and then the corresponding mass- 
effective temperature {Tcs) diagram in Figure [TTJ 

We also perform a qu antitative com parison with the 
same models. We follow pTorresI (|2007t ) in defining a pa- 
rameter equal to the ratio of the observed radius to that 
predicted by the models given the observed mass, al- 
though we denote this parameter by p instead, to avoid 
any potential confusion with the (3 parameter defined by 
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TABLE 11 

Detached, double-lined, double-eclipsing main sequence EB components below 0.4 Mq^ 



Name 


Period 


M 


R 


(K) 


[M/H] 


Source'' 




(days) 


(M©) 


(Rq) 




oiJoo-iViiljJD-1 A 


U.4U ( 


n 070 J- n non 
yJ.Z 1 Z ± U.UzU 


U.ZDo it U.UUyU 






1 

i. 


OiJoo-iVilljJD-i ID 




n 0/in J- n noo 
U.z4U it U.yJZZ 


U.z4o it U.UUo4 


ooUU it loU 






T QOQ^i A C 
LrJ OZOO A 


n v'zi 
a. 1 1 L 


U.o/D it U.Uiu 


U.vSiyo it U.UUo4 


1 o 1 1 1 n 
ooLZ ± liU 




o 
z 


P T '^9'?*^ 
\jo OZOU -D 




n OQl -1- 01 


U.OUU ^ U.U-LO 








CM Dra A 


1.27 


0.2310 ±0.0009 


0.2534 ±0.0019 


3130 ± 70 


[-1,-0.6] 


3 


CM Dra B 




0.2141 ±0.0010 


0.2396 ± 0.0015 


3120 ± 70 




LP 133-373 A 


1.63 


0.340 ±0.014 


0.33 ±0.02 


3058 ± 195 




4 


LP 133-373 B 




same 


same 


3144 ± 206 






MGl-2056316 B 


1.72 


0.382 ± 0.001 


0.374 ±0.002 ±0.002 


3320 ± 180 




5 


KOI 126 B<1 


1.77 


0.2413 ±0.0030 


0.2543 ± 0.0014 




+0.15 ±0.08 


6 


KOI 126 




0.2127 ±0.0026 


0.2318 ±0.0013 








CCDM J04404-f3127 B.C^ 


2.05 










7 


CU Cnc B 


2.77 


0.3980 ±0.0014 


0.3908 ± 0.0094 


3125 ± 150 




8 


IRXS J154727.5-I-450803 A 


3.55 


0.2576 ±0.0085 


0.2895 ± 0.0068 






9 


IRXS J154727.5-I-450803 B 




0.2585 ±0.0080 


same 









^ This mass criterion has been apphed to keep the number of objects in the table and plots manageable, 
and is not intended to necessarily be physically meaningful beyond being appropriate for comparison 
to the present s yste r n. 

b (1 ) 31ake et ap ([gOOSl ), (2) llrwin et al.l (l2009bD. (3) IMorales et aH fl^Oj (4) IVaccaro et aP 
(pool. (5) iK raus et al. (20l l. (61 [Carter et al\mm . (7) IShkolnik et al.l (120101) . (8) iRibasI (|20Q1 : 
iDelfosse et al., (,199a) . f9) .Hartman et al.l (|20lH ). 

° Parameters determined giving equal weight to all three models, following iHartman et al.l (|2011| ). 
<i While not double-lined, this object is a special case as it is still possible to solve for the masses and 
radii of both M-dwarfs independent of IVI-dwarf models. The period given is that for the inner M-dwarf 
binary as this is presumably the appropriate one for estimation of activity levels. 
<= Full solution not available. 



iChabrier etUI (|2007f ) and used bv Morales et al.l (|2010D . 
which is not the s ame. 

As argued in Torres! ()2007f) . there is evidence that the 
bolometric luminosities from the models are not seri- 
ously in error. If this is the case, a corresponding fac- 
tor of p^^/^ must be applied to the effective tempera- 
tures from the models when inflating the radii, in order 
to preserve the bolometric luminosity. We show results 
for this scenario in the tables and plots. Note, how- 
ever, that the bolometric lumino sity is not preserved 
in several recent theoretical works (IChabrier et al.|[2007l : 
IMorales et aI1l2010l: MacDonald fc Mullanll201lD . so the 
temperatures derived from this procedure should not be 
considered definitive. 

As is often the case in eclipsing binaries, while the indi- 
vidual masses of the objects are well-determined from the 
radial velocities, the sum of the radii is much better de- 
termined than the individual radii, modulo limb darken- 
ing uncertainties as discussed in !j4.5l We therefore com- 
pute, in addition to the p values for the individual stars, a 
"composite" value comparing the model-predicted radius 
sum given the individual observed masses, to the mea- 
sured radius sum. This quantity, which we call pi+2 , has 
much smaller formal uncertainties than the individual pi 
and p2 values, and measures essentially a "weighted av- 
erage inflation" for the whole system. These quantities, 
and the corresponding effective temperatures predicted 
from preserving the bolometric luminosity, are given for 
a range of assumptions for the age and metallicity in Ta- 
ble [121 While we caution over- interpretation of the pi+2 
results due to limb darkening uncertainty, the individual 
values should be much less prone to this problem. 



5.1. Mass-radius comparison 

We first limit our comparisons to the mass-radius 
plane. If the metallicity is solar, the sum of radii is in- 
flated at the 7(7 level compared to the models, even for 
10 Gyr age, which produces the largest model radii. The 
discrepancies are not as significant in the individual radii, 
with the primary 2a and the secondary Scr larger than 
the model predictions. We also note that the "eclipsed 
spots" model on the secondary, which we discarded in 
M.'Sl does not solve the radius discrepancy, having a sec- 
ondary radius (and sum of radii) even further from the 
model predictions. 

Youth is unlikely as an explanation given the old disk 
kinematics of our target. We estimate an age of approx- 
imately 120 Myi at solar metallicity would be needed 
to bring the radius sum into agreement with the models, 
yet such objects should also rotate extre mely rapid l y and 
exhibit strong Ha em ission (e.g. Terndrup et al.|[2000l: 
IHartman et al.l [20101 and references therein), neither of 
which are seen. 

Metallicity may provide a viable explanation for the 
observed radii given the present uncertainties (both in 
the observations and the models, where for the latter it 
is possible the effect of metallicity on ra dius is under- 
predicted; see ilT] and [Lopez- Moralesll20(57l ) . although the 
object would most likely need to be extremely metal rich. 
This would place it in the tail of the metallicity distribu- 
tion for the old Galactic disk populat ion, which is li kely 
to peak around [M/H] w -0.5 (e.g. iLeiggettl [19921 and 
references therein). 

Given that metallicity has difficulty reproducing radii 
in other well-observed systems (see Sjl} a-nd the kinematic 
evidence, we proceed to examine other possible causes of 
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TABLE 12 

Quantitative mass-radius comparison of LSPM J1112+7626 with models. 



Age 


[M/H] 


Pi 


P2 




Tea. 2 


T'cff,2/T'cff,l 


(Gyr) 




Pl+2 


(K) 


(K) 




1 


0.0 


1.054 ±0.016 


1.096 ±0.018 


3426 ± 27 


3255 ± 26 


0.950 ±0.013 


1 


-0.5 


1.070 ±0.016 


1.130 ±0.018 


3621 ± 28 


3387 ± 27 


0.935 ± 0.013 


5 


0.0 


1.037 ±0.015 


1.072 ±0.017 


3453 ± 27 


3293 ± 27 


0.954 ±0.013 


5 


-0.5 


1.055 ±0.016 


1.101 ±0.018 


3648 ± 29 


3436 ± 28 


0.942 ± 0.013 


10 


0.0 


1.029 ±0.015 


1.052 ±0.017 


3468 ± 27 


3324 ± 27 


0.959 ± 0.013 


10 


-0.5 


1.047 ±0.016 


1.078 ±0.017 


3666 ± 29 


3475 ± 28 


0.948 ±0.013 



1 0.0 l-0™6lQ;oo5g 3396 ± 13 3291 ± 12 0.96906 ± 0.00035 

1 -0.5 l-094llg;°°|^ 3579 ± 14 3439 ± 13 0.96096 ± 0.00045 

5 0.0 l-0505to.oo53 3428 ± 13 3323 ± 12 0.96940 ± 0.00038 

5 -0.5 l-0735lo;o°55 3614 ± 14 3477 ± 13 0.96195 ± 0.00050 

10 0.0 1.0378lg;Q°|2 3451 ± 13 3344 ± 12 0.96915 ± 0.00040 

10 -0.5 l-0590to °°56 3642 ± 15 3504 ± 13 0.96221 ± 0.00052 



the inflated radii. 

It is tempting to suggest that the inflation in 
LSPM J1112+7626 results from elevated activity levels 
in one or both components. While the lack of Ha emis- 
sion, unless this is rendered undetectable by the limited 
signal-to-noise ratio of our spectra, may be a difficulty 
for this scenario to explain, it should be noted that more 
modest levels of activity can result in weakened Ha emis- 
sion, or even absorption, but may still produce sufficient 
spot coverage to explain the inflation. 

This system straddles (in mass) the full convection 
limit, where a substantial increase in t he observed activ - 
ity lifetimes for M-dwarfs occurs (e.g. IWest et al.l[2008l) . 
This would argue that the secondary is more likely to 
be active and thus the source of the inflation. The sec- 
ondary lines in the spectrum are also more difficult to 
observe due to its lower luminosity, which would make 
an Ha line easier to hide in the noise. Identification of 
the star hosting the out of eclipse modulations may shed 
additional light on which of the components might be 
responsible for the inflation. 

As discussed by Morales et al.l (|201Cll ). the radius in- 
flation inferred from eclipsing binary analyses often im- 
plies extremely high spot filling factors. We estimate that 
p — 1.05 corresponds to /3 = 0.1. Assuming a spot tem- 
perature contrast of Tcff^spot/Tcfi^phot = 0.9, this would 
require a filling factor of approximately 30%. This is 
much larger than needed to produce the observed out 
of eclipse modulation, and would mean there is a sub- 
stantial longitudinally symmetric spot component, if the 
inflation is indeed caused by spots. We showed in M.3.31 
that such a level of spots does slightly modify the com- 
ponent radii if they are not eclipsed, reducing R2 and 
increasing Ri. In practice, this may slightly reduce the 
required filling factor depending on where the spots are 
located, but unless the spot coverage is much larger than 
we assumed, the systematic errors in the radii resulting 
from such spots are still insufficient to explain the obser- 
vations without invoking inflation. 

Finally, on suggestion of the referee, we note that this 
object may indicate the need to revisit the equation of 
state for very low mass stars, which was discussed as a 
possibility to explain the radius inflation before the ac- 
tivity hypothesis gained prevalence (e.g. .Lopez-MoraleSi 



120041: [Torres fc Ribasll2002[ ). The equation of state was 
a significant source of difficulty in early work attempt- 
ing t o model very low mass stars (Chabricr & BarafQ 
119951 ). and although substantial progre ss on this front 
has been made (see the r eviews by iChabrier fc Bara!^ 
I2OOOI: IChabrier et al.ll2005[ ). it is still an open question. 

5.2. Mass-Tcfi comparison 

Comparisons in only the mass-radius plane ignore im- 
portant information contained in the effective tempera- 
tures. The correct physical model must explain all of the 
observations simultaneously, so we now proceed to exam- 
ine the temperatures. This comparison is more problem- 
atic than in mass-radius for a number of reasons. The 
main observational issues are the difficulty of constrain- 
ing effective temperatures (see ^14. 81) and metallicities for 
the eclipsing binaries. 

As noted in Sjl] unlike the radius, the effective temper- 
ature predicted from models does depend quite strongly 
on metallicity, because the bolometric luminosity is a 
function of metallicity. As the metallicity is decreased, 
the model bolometric luminosity and effective temper- 
ature increase, and the radius decreases by a small 
amount. While metallicity does complicate the inter- 
pretation of the mass-Toff diagram, this a key argument 
against metallicity as the explanation for all of the in- 
flated radii in the eclipsing binary sample, as discussed 
already in SjTl In. this regard, CM Dra is puzzling as 
several authors have claimed this object is metal poor, 
which would further exacerbate the effective temperature 
discrepancy with the models. 

Given the present uncertainties, it is not clear 
if the effective temperature difference between 
LSPM J11124-7626 and the models is signiflcant. 
However, it does lie at the low-temperature end of the 
eclipsing binary results, and this lends support to the 
hypothesis that this object is more metal rich than the 
other eclipsing binaries. 

5.3. Future work 

It is clear that LSPM J1112+7626 offers one of the best 
prospects to observationally test the causes of inflated 
radii in eclipsing binaries, but further observations are 
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needed. The most fruitful avenues would be to pursue 
a determination of the metallicity, improvement of the 
system parameters with a particular focus to investigat- 
ing potential systematics in the radii (e.g. due to spots), 
and constraining the activity levels in the components by 
independent means (for example, X-ray emission, or the 
Can H and K lines). 

While we have mentioned the non-synchronized spin 
and orbit, we were unable to take full advantage of this 
property with the available observational data. As shown 
in the effect of spots should be larger on the primary 
eclipse, so we advocate an intensive campaign on these 
events, preferably to obtain multiple, complete eclipses 
with the same telescope and detector system at a range 
of rotational phases and in multiple wavelengths. The 
primary eclipses are observable from Europe and Scan- 
dinavia during the winter 2011-2012 season, where ex- 
cellent observational facilities are available. We remind 
observers that good-quality out of eclipse monitoring will 
be important for the interpretation of the eclipses them- 
selves, and the spot distributions are likely to evolve (in- 
deed. Figure [7] indicates that we may have already seen 
this in the existing data), so the monitoring must be 
contemporaneous, and preferably also taken at multiple 
wavelengths to constrain the spot temperatures. 
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